Phase transitions and equilibrium measures in random matrix 

models 



A. Martmez-Finkelshtein, R. Orive and E. A. Rakhmanov 
^ March 5, 2013 

o 

Abstract 

. The paper is devoted to a study of phase transitions in the Hermitian random matrix 

models with a polynomial potential. In an alternative equivalent language, we study 
families of equilibrium measures on the real line in a polynomial external field. The total 
mass of the measure is considered as the main parameter, which may be interpreted also 
either as temperature or time. Our main tools are differentiation formulas with respect 
to the parameters of the problem, and a representation of the equilibrium potential in 
terms of a hyperelliptic integral. Using this combination we introduce and investigate a 
dynamical system (system of ODE's) describing the evolution of families of equilibrium 
measures. On this basis we are able to systematically derive a number of new results on 
phase transitions, such as the local behavior of the system at all kinds of phase transitions, 
as well as to review a number of known ones. 
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1 Introduction 

This paper is devoted to a study of families of equilibrium measures on the real line in 
the polynomial external fields. We consider these measures as functions of a parameter t 
representing either the total mass, time or temperature. These families are regarded as 
models of many physical processes, which motivates their intense study in the context of 
the mathematical physics. But equilibrium problems for the logarithmic potential play an 
important role also in analysis and approximation theory; in particular, they provide a general 
method used in the theory of orthogonal polynomials. In this sense, the subject is essentially 
"bilingual" and has so many ramifications that we opted for writing an extended introduction, 
instead of a dull enumeration of our results. 

One of our main goals is a formal investigation of the dependence of the key parameters 
of the equilibrium measure from the parameter t, and in particular, their singularities as 
functions of t. In thermodynamical terms, these singularities are closely related to phase 
transitions, and we start the introduction with the discussion of this problem in the context 
of statistical mechanics and thermodynamics. It does not mean that we actually intend to 
interpret our results in any nonstandard way. We want, on the contrary, to make a short 
review of a few existing standard interpretations, which could be helpful during "translations" 
in our bilingual area. 

Then we formally introduce basic facts related to the equilibrium measure and mention 
briefly some applications in analysis and approximation theory, returning at the end to the 
topic of the Coulomb gas and random matrices. 

In the course of these discussions we describe in general terms the main results of the 
paper and its structure. 

1.1 Statistical mechanics 

Determinantal random point processes (or self-avoiding point fields) are pervasive in sta- 
tistical mechanics, probability theory, combinatorics, and many other branches of mathe- 
matics [161 EH US S3 E3 E2] • Among the examples of such processes we can mention 2-D 
Fermionic systems or Coulomb gases, Plancherel measures on partitions, two-dimensional 
random growth models, random non-intersecting paths, totally-asymetric exclusion processes 
(TASEP), quantum Hall models, and random matrix models, to mention a few. Many of these 
models lead to the so-called orthogonal polynomial ensembles [48, 82j, an important subclass 
of determinantal random point processes. The large scale behavior of such ensembles is de- 
scribed in terms of the asymptotics of the underlying family of orthogonal polynomials, and 
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in the last instance, by the related equilibrium measure solving an extremal problem from 
the potential theory. 

One of the earliest and probably best known examples of orthogonal polynomial ensembles 
is the joint distribution of the eigenvalues of a random N x N Hermitian matrix drawn from 
a unitary ensemble (see [32] ) . More precisely, we endow the set of N x N Hermitian matrices 

[M = (M jk )» k=1 : M kj = M~k} 

with the joint probability distribution 

^ N N 

dv N (M) = ^- exp (-Tr V(M)) dM, dM = JJ dMjj \\dRe M jk dlm M jk , 



j=i j^ k 



where V : K — )■ R is a given function with enough increase at ±00 to guarantee the conver- 
gence of the integral in the definition of the normalizing constant 

Zn = j exp(-TrV(M)) dM. 

Then vn induces a (joint) probability distribution \in on the eigenvalues Ai < • • • < Xn of 
these matrices, with the density 

7 ^ A ) = Yn~ - A i) 2 ex p ( - E . W 

i<j \ i=l / 



Vn 

and with the corresponding partition function 
Zn 



I ... I WiX, - Xj) 2 exp f - V V(Xi) \ d\ 1 . . . dX N . 



i<3 

This result can be traced back to a classical theorem by H. Weyl, see |83j and also [60] . Notice 
that vn is invariant under unitary transformations of the Hermitian matrices, while measure 
jjiN defines the equilibrium statistical mechanics of a Dyson gas of N particles with positions 
Xj on M in the potential V, where they interact by the repulsive electrostatic potential of the 
plane (the logarithmic potential). On the other hand, is also the joint probability density 
function of an iV-point orthogonal polynomial ensemble with weight w(x) = exp(— V(x)) (see 
e.g. \28\ EE]). The free energy of this matrix model is defined as 

Fn = ~^logZ N - 

Regardless its interpretation, a particularly important problem is to analyze its asymptotic 
behavior in the thermodynamic limit, i.e. as N — > 00. 

A rather straightforward fact, the existence of the limit 

Foo = lim F N 

N^oo 
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(infinite volume free energy) has been established under very general conditions on V, see 



e.g. [16] and Section 1.2 below. An important property of the infinite volume free energy 
is its analyticity with respect to the parameters of the problem. The values of the parameters 
at which the free energy is not analytic are the critical points; curves of discontinuity of some 
derivatives of the free energy connecting the critical points divide the parameter space into 
different phases of the model. Thus, critical points are points of phase transition, and they 
are going to be a center of our case study. 

The observation that the distribution ([I]) can be regarded as the Gibbs ensemble on the 
Weyl chamber {A: Ai < • • • < Ajy} with Hamiltonian 

1 1 N 

£ log K^ + 2^>^ (2) 

i<j ISl «i j=i 

allows to foretell the fundamental fact that the value of -Fqo is given by the solution of a 
minimization problem for the weighted logarithmic energy. The corresponding minimizer is 
the equilibrium measure associated to the problem^ This measure, which is a one-dimensional 
distribution on E, is also a model for the limit distribution of the eigenvalues Ai < • • • < Ajv 
in 0. Indeed, the multidimensional probability distribution fi]\f is concentrated for large 
N near a single point A* = (A*,...A^) 6 Mr, which is actually the minimizer for the 
corresponding discrete energy Q of an Appoint distribution. In other words, for large N the 
measure /ijy is close to <5(A — A*). As N — > oo, the discrete equilibrium measure converges 
to the continuous one, so that the thermodynamic limits are essentially described by the 



(continuous) equilibrium measures (see Section 1.3). 

When V is real-analytic, the support of such a measure (or the asymptotic spectrum 
of the corresponding unitary ensemble) is comprised of a finite number of disjoint intervals 
(or "cuts"), and the number of these intervals is a fundamental parameter. For instance, 
the free energy i*jy has a full asymptotic expansion in powers of N~ 2 ( "topological large N 
expansion") if and only if this support is a single interval; otherwise, oscillatory terms are 
present (this was observed in [15] , and studied systematically in [131 E21 EH] ) . 

The particularly interesting phenomena occur precisely in the neighborhood of the values 
of the parameters at which the number of the connected components of the support of the 
equilibrium measure changes. Any change in the number of cuts is a phase transition in the 
sense specified above, but there are also phase transitions of other kinds (not related to a 



change in the number of cuts): see Sections 3.3 and [4] for details 



We specialize our analysis to the polynomial potential. It is convenient to write V in the 
form 

V(x) = V n (x) = 2n<p(x), (3) 

where (p is a polynomial of an even degree and positive leading coefficient. This case is of 
great interest, see e.g. [H E] HU EI] , to cite a few references. A common situation is when 
n, N — > oo in such a way that 

N 

lim — = t > 0. (4) 

n 



x We have gathered the related definitions and basic facts in Section 1.2 



4 



Although the coefficients of <p ( "coupling constants" of the model) play a role of the variables 
in the problem, the parameter t > stands clearly out, as it was already mentioned above. 
Recall that it can be regarded either as a temperature (from the point of view of statistical 
mechanics) or time (from the perspective of a dynamical system), and will correspond to the 
total mass of the equilibrium measure in the external field ip on M. 

One of the goals of this paper is the description of the evolution of the limiting spectrum 
of the unitary ensemble for Q-Q as time (temperature) t grows from zero to infinity, paying 
special attention to the mechanisms underlying the increase ( "birth of a cut" ) and decrease 
("fusion of two cuts" or "closure of a gap") in the number of its connected components. 
Equivalently, we study one-parametric families of equilibrium measures on R in a polynomial 
external field with the total mass t of the measure as the parameter. 

This problem has been addressed in several publications before, here we only mention 
a few of them. For instance, [9] studied the fusion of two cuts for the quartic potential 
<p, proving the phase transition of the third order for the free energy (in other words, that 
-Foo G C 2 \ C 3 ). This work found continuation in [U[36], where in particular the birth of the 
cut in the same model was observed (but not proved rigorously). It turns out actually that 
the third order phase transition of the infinite volume free energy is inherent to all possible 
transitions existing in this model. 

The evolution of the support of the equilibrium measure in terms of the parameters of 
the polynomial potential has been studied also in [3], SI 132] , where a connection with some 
PDE's and integrable systems have been exploited. 

However, despite of this intense activity, the picture is not complete; even the monographic 
chapter [22J contains imprecisions. Open questions exist actually for the first non-trivial case 
of a quartic field, and this paper adds a number of new details to this particular picture; 
some of them seem to be significant. For instance, we present a simple characterization of 
the case when the equilibrium measure has one cut for all values of the parameter t > 0, or 
when a singularity of type III (see the definition in Section 1.4) can occur. We also study in 
more detail the system of differential equations governing the dynamics of the endpoints of 
the support of the equilibrium measure. In particular, we analyze the behavior of the infinite 
volume free energy and of other magnitudes of the system near the singular points, revealing 
some interesting universal properties. 

Finally, but not less important, we show that all the known and new results related to the 
outlined problem can be systematically derived from two basic facts in the potential theory. 
One of them is a representation for the Cauchy transform of the equilibrium measure. This 
representation is a corollary of the fact that on the real line any equilibrium measure in an 
analytic field is a critical measure, which allows us to apply a variational technique presented 
in [3D], and systematically in [58] . Another fact is a Buyarov-Rakhmanov differentiation 
formula [21] for the equilibrium measure with respect to its total mass and some of its 
immediate consequences, which we complement with a unified treatment of the differentiation 
formulas with respect to any coupling constant, a result which to our knowledge is new. We 
provide further details in Section [pj 



5 



1.2 Equilibrium measure in an analytic external field on R 

In this section we introduce basic notation and mention a number of fundamental facts on 
the equilibrium measures on the real line, necessary for the rest of the exposition. For more 
details the reader can consult the original papers [211 E3 EH EH [691 EO] arid the monograph 



For a finite Borel measure a with compact support supp(cr) on the plane we can define 
its logarithmic potential 

V cr (z) = — J log \z — x\ da(x), 

and its logarithmic energy, 



I [a] = — J J log \z — x\ do{x)do(z). 



Suppose further that a real- valued function ip, called the external field, is defined on supp(cr). 
Then we introduce the total (or "chemical") potential, 

W°(z) = V%z) + p{z), (5) 

(defined at least where ip is) and the total energy, 

I lp [a]=I[a}+2 J <p(z)dz, (6) 

respectively. This definition makes sense for a very wide class of functions <p, although for the 
purpose of this paper it is sufficient to consider basically real-analytic, actually polynomial, 
external fields on the real axis. As usual, we assume also that 

lim (p(x)/log \x\ = +oo, 

\x\— >+oo 

condition that is automatically satisfied for any real non-constant polynomial of even degree 
and positive leading coefficient. Then for each t > there exists a unique measure At = 
Xt((f,M) with compact suport St = supp(At), minimizing the total energy 

I v [Xt] = mm I v [<r] (7) 

crGMt 

in the class Ait of positive Borel measures a compactly supported on M and with total mass t 
(that is, o~(M) = t). Moreover, Xt G Ait is completely determined by the equilibrium condition 
satisfied by the total potential 

Wi>(z)l =Ct > ^ SUPP ^' (8 ) 
v \> ct, zeM, 

where the equilibrium or extremal constant a can be written as 

d = ct(<p) := j Up [At] - f ^d\A • (9) 
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When K is a compact set on M and the external field cp is 



<p(x) 



0, if x G if, 
+00, otherwise, 



the corresponding energy minimizer Xt(ip,M) in the class of all probability measures A4\, 
denoted by ojk, is known as the Robin measure of K, its energy is the Robin constant of K, 

p(K) = I[oj K ], (10) 

and the logarithmic capacity of K is given by cap(K) = exp(— p{K)). Measure lok is char- 
acterized by the the fact that supp(w/<-) = K together with the equilibrium condition 

V UK (z)=p(K), zeK. (11) 

For a general external field the main technical problem when finding the equilibrium 
measure is that its support is not known a priori and has to be found from the equilibrium 
conditions. Once the support is determined, the measure itself can be recovered from the 
equations presented by the equality part in ([8]), which is an integral equation with the log- 
arithmic kernel. After differentiation it is reduced to a singular integral equation with a 
Cauchy kernel whose solution (on reasonable sets) has an explicit representation. 

Hence, solving the support problem is the key, and it is essentially more difficult. 

In the one cut case, its endpoints are determined by relatively simple equations, but in 
general, straightforward solutions are not available. For a polynomial (actually, real analytic) 
field the support St is a union of a finite number of intervals. Equations may be written for 
the endpoints of those intervals (and we better know the number of intervals in advance), 
but these equations are not easy to deal with. For instance, in the case of the polynomial 
field they may be interpreted as systems of equations on periods of an Abelian differential on 



a hyperelliptic Riemann surface (see e.g. Remark 3.2), but such systems are usually far from 
being simple. In this context, our approach is based on the dynamics of the family St- As we 
have mentioned above, one of our goals is to describe the dynamics of the family of supports 
St as the mass t (which is also "time" , or according to [9] , the "temperature" ) changes from 
to +00. The detailed discussion starts in Sections [2] and [3] below. 

Observe finally that a study of the family with total mass t as a parameter in a 
fixed external field 92 may be reduced by a simple connecting formula to the family of unit 
equilibrium measures with respect to family of fields jtp. We have 

c t (<p)=tc 1 (j<p\ \ t (v)=t\i(^ip). (12) 

In some cases such a reduction is reasonable, but more often it brings more problems than 
benefits. In particular, we believe that the total mass as a parameter has particular advan- 
tages in the problem under consideration. 
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1.3 Weighted equilibrium in Analysis 

The significant progress in the theory of rational approximations of analytic functions (Pade- 
type approximants) and in the related theory of orthogonal polynomials in the 1980's is due 
in part to the logarithmic potential method; see the original papers [381 EH1 EH EH |65j E3 
EHl H3 EH E3 EE] and the monographs [EH E2J E2], although the list is not complete. A 
novel important ingredient of the techniques developed during that period was precisely the 
introduction of the weighted equilibrium measure (among other types of equilibria). 

For instance, a typical application of the notion of equilibrium measures for orthogonal 
polynomials is the following fundamental result from [39J (see also |544 loT] 172]). presented 
here in a simplified form: let ip be a continuous external field, re, N S N, and let polynomial 
Qn{x) = Qn,u{x) = x N + . . . be defined by 

Q N {x)x k e- 2n ^ x) dx = 0, k = 0, 1, . . . , N - 1. 

This polynomial is also characterized by the extremal property 

m N : = f Q 2 N {x)e- 2n ^ x) dx = min [ P 2 {x)e~ 2nip{ - x) dx. 

JR P{x)=x n +... J R 



Now, if n, N — > oo in such a way that N/n — > t > 0, then 

±lo g (m N )=c t , - 



lim lo s( m ^v) = c *> Jj Yl ^C^Hf), (13) 

Qiv(C)=o 



where — ->■ is the weak-* convergence. From the point of view of the spectral theory, zeros of 
Qn are eigenvalues of the truncated N x N Jacobi matrices associated with the weight e~ 2n<fi , 



and (13) shows that Xt(f) can be naturally interpreted as the limit spectrum of the infinite 
Jacobi matrix, associated with the scaling N/n — > t. For details related to applications of 
equilibrium measures in spectral theory of discrete Sturm-Liouville operators see |75j . 

Another large circle of applications of (continuous) equilibrium measures Xt(f) is related 
to the discrete analogue of this notion. Let 

M n := |m = 5^C*' CfcSlR, fc = l,...,nj C M n 

denote the set of all point mass measures on K of total mass re. The corresponding discrete 
energies are defined by 



£[/i]:=5>g— ^— , := E\pt] + 2 f 



ip dfj, 



(cf. (§). 

Let [i n = finif) ^ M n be a minimizer (not necessarily unique) for E^lfj] in the class M n \ 
points {Ci,n> • • • , Cn,n} = supp(/i„) are also known as weighted Fekete points. We have 



1 * 

- n n (ip) — > Ai(^), 
re 



1 

n 



as n — > oo; 



S 



for the proof, see e.g. [72J, where some applications of Fekete points in approximation theory 
are also discussed. Furthermore, for a class of external fields cp, the polynomials 

n 

H n {x) Yi ( T ^ fc > n ) 
fc=i 

satisfy a second order differential equation with polynomial coefficients and constitute an 
example of the so-called Heine-Stieltjes polynomials; see (56J E3 EE1 EHJ E7J E2 El] for further 
details on the algebraic and analytic properties of these polynomials. As it was shown in [59], 
the limit zero distribution of such polynomials in the classical case studied by Stieltjes [50] 
is governed again by the (continuous) one-parametric family of equilibrium measures. Also 
a related work worth mentioning is [B], where deformation with respect to the parameters of 
the external field were considered. 



1.4 Random matrices and Dyson gas revisited 

Let us return to the random models mentioned in the beginning of this Introduction. An 
important fact, already mentioned there, is that the eigenvalue distribution of the unitary 
ensemble of random matrices ([!]) may be equivalently written in terms of the discrete energy. 
Indeed, since we deal only with symmetric functions of a vector £ = (Ch • • • > Cn) G we 
may identify £ with the non-ordered collection { (j,. . . ,£ n } of real points. This collection 
uniquely determines measure rj 6 M n with supp(^) = {Ci, • • • , Cn} by 

n 

i=l 

so that we can identify also £ £ W 1 rj £ M n and write E v (£) instead of E^rj). 
Now, we define 

dVn,N(0 = ^-e- n2E ^\ Z niN = [ e~ n2E ^dC (14) 

where if) = ^<P- Distribution jU n ,jv(C) coincides with cI[j,n(\) in ([T]) if we identify A = £ 
and use also connection ([3]). As observed before, (14) is a canonical Gibbs distribution for a 
Dyson gas. 

In this connection, the classical Heine's formula 

N 



Qn(x)= / TT(x - Cfc) djl n ,n(C), 

Jr n , , 



=i 

represents the orthogonal polynomial Qn with respect to the weight e~ Nip = e~ 2n ^ as the 
average of polynomials with random real zeros distributed according to /x ni jv- 

In the last decade the distribution n n ,N and its characteristics, such as the correlation 
function, has been intensively investigated in relation to the eigenvalue statistics of the unitary 
ensembles of Hermitian matrices. 
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One of the relevant features of /J, n N i n the asymptotic (thermodynamic) regime is that 
it exhibits non-trivial phase transitions or critical behavior, with universal properties that 
attract increasing attention. The existence of such phase transitions and their character has 
important physical interpretations, see e.g. [2"1 129] . 

The equilibrium measure At (and consequently, the parameters upon which it depends) is 
called regular (see e.g. |27j ) if the following conditions hold: (a) Af is absolutely continuous 
on R with respect to the Lebesgue measure and X' t does not vanish in the interior of St', (b) 
X' t has a square root vanishing at the endpoints of the connected components of the support 
Sf, and (c) the inequality in ^ is strict in K \ St- Otherwise, the equilibrium measure is 
singular. Following the terminology of \27\ I51j . an equality in ^ at a point in K \ St is a 
singularity of type I (this is also known as "birth of the cut" , see below) ; zeros of X' t that lie 
in the interior of St are called singular points of type II (corresponding to a closure of a gap), 
while the endpoints of the connected components of St where X' t has a higher order vanishing 
are called singular points of type III. 

It is rather easy to see that the infinite volume free energy or the correlation kernels are 
analytic functions of the parameters of the problem as long as the underlying equilibrium 
measure is regular. Hence, phase transitions (in the physical sense of break of analyticity 
of the free energy) are intimately connected to the singularities of the equilibrium measure; 
this is the central idea of this work. Our goal is to discuss the mechanisms of these phase 
transitions and to analyze the behavior of the free energy near its critical values. 

The asymptotic behavior of a random matrix ensemble near singular points of the equi- 
librium measure has been studied intensively (see [TUJ, [HI [231 [231 ED EH S3 [62] , to mention 
a few). A remarkable feature of this behavior is revealed in the so-called double scaling limit, 
where the parameters near the critical value are coupled with the N (the size of the matrix) . 
In this case, with an appropriate scaling of the variables, the asymptotic behavior of the cor- 
relation kernel near the critical point depends only on the type of the singularity rather than 
on the potential itself, which is one of the manifestations of the universality in the random 
matrix models. 

It is known that the asymptotic behavior of the correlation kernel in a double scaling 
limit at the singular points is connected to the hierarchies of the Painleve equations. In the 
general polynomial case, the local behavior at the birth of the cut was described in [33] (in 
terms of orthogonal polynomials with respect to the Freud weight e~ x " on K), and proved 
almost simultaneously in [71 [23"1 [62] using the Riemann-Hilbert techniques. For the type II 
singularity [111 [2"%] , the kernel can be written in terms of the Hastings-McLeod solution of the 
Painleve II equation, while for the type III singularity, the kernel is related to the Painleve I 
transcendent [23 EI] • The double-scaling asymptotic analysis of this sort is beyond the scope 
and techniques of this paper. 

2 Families of the equilibrium and critical measures on the real 
line 

In the study of measures minimizing the total energy functional ^ it is often convenient 
to consider simultaneously solutions of certain more general but closely related types of 
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equilibrium problems and associated measures, which we call critical measures. We discuss 
here two kinds of critical measures - the local minima and the saddle points of the energy, 
used to derive an equation for the Cauchy transform of equilibrium measure (see Theorem 



2.2 below). In this connection we make also a few general remarks on the equations that can 



be useful for the constructive determination of the equilibrium measures. 
2.1 Integral equations for the density 

We fix a smooth external field tp G C 1+£ (IR) and t > 0. Recall that the pair consisting of a 
positive measure A = At G A^t(K) (the equilibrium measure with total mass t) and and an 
(equilibrium) constant c = q £ 1 is uniquely defined by the (equilibrium) conditions Q; 
measure At is also the unique global minimizer of the total energy in the class .Mi(ffi). Under 
the assumptions above, it is absolutely continuous with respect to the Lebesgue measure, so 
that dX(y) = \'{y)dy on R. 

Consequently, the problem of the constructive determination of At boils down to finding 
its support St, a problem that has been under investigation for a few decades. The main 
conclusion is that this problem has a simple constructive solution if St is a single interval 
(see e.g. [53] or [721 Ch. IV]); when the support is not connected, a satisfactory constructive 
solution is not available. 

Let us provide some details in this respect. Assume that the support is a union of p > 1 
disjoint intervals, 

p 

St = (J A*,, A fc = [a2k-i,a>2k], a 1 < a 2 < ■ ■ ■ < a 2p - (15) 
fc=i 

Differentiating the equality part of the equilibrium conditions ^ one obtains the singular 
integral equation 

A'(y) 



s t x ~y 



dy = (p'(x), x G S t . (16) 



This equation on the unknown function A' G L 1 (St), given a fixed St in (15) and a fixed 
if' G C £ (St), is a well known object in the theory of singular integral equations. A general 
solution of this equation may be written in terms of singular integrals (see [H Ch. 5, §4], 
as well as \37\ I63j). The explicit expression for this solution is not so relevant here, but it 



renders a system of p equations on the 2p endpoints aj of the support St (see equations (17) 
below) that are central to our discussion. 

Condition ip G C 1+e (R) implies that A^ G C(St), and in particular, that the density of the 
equilibrium measure is bounded on St- Thus, we are interested in the bounded solutions of 



(16), and by [8], it yields the following result: 



Theorem 2.1. Assume t > fixed. A bounded solution for (16) exists if and only if 

j '( \ 2p 
1J 7 == dy = 0, j = 0, . . . ,p - 1, where A(z) = JJ (z - a k ). (17) 
St y\A{y)\ k=l 
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Furthermore, in this case the bounded solution is unique and is given by the formula 



x'iv) = -v\m\ 

IT 



dx. 



(18) 



I St (x-y)y/\A(x)\ 

As we have mentioned above, this theorem is a corollary of a representation for general 



solution of (16) [8] |371 E3] ; an alternative derivation of equations (17) may be found in 
Ch. 6]. 



A general solution of (16) satisfies 



V Xt (x) + ip{x) 



-t > 



x G At, k = 1, 2, 



,p, 



(19) 



(k) 

where constants c t are not necessarily the same on different intervals Afc. In order to get 
(Tsl) we must impose p — 1 extra conditions, 



,(1) 



.(2) 



» 



(20) 



Since Aj is uniquely determined by tp and St, these are equations on a^s. Actually, equations 
(17) and (20), complemented with the normalization condition \t(S) = t, constitute a system 
of 2p equations on 2p unknowns, the endpoints of the intervals in St- Using the expression 



for At in (18) they can be explicitly rewritten in terms of ip. We omit these calculations here, 



but see (24) below, where these equations will be rewritten in a form more suitable for our 
needs. 

The question in what sense these equations determine St = supp(At) is not trivial. Actu- 
ally, we still need to verify the inequalities X'(y) > 0, y G St, along with V x (x) + tp(x) > ct, 
x € M. If they are satisfied, the problem is solved. Can this solution be regarded as con- 
structive? We can certainly assert this in the one-cut situation. But we are inclined to give 
a rather negative answer in the multi-cut case. To the best of our knowledge, the general 
support problem is actually a whole research program. 

Let us point out that once St has been established, it determines the partition of the total 
mass t into its "portions" 



I At | (A, 



d\ t , k = l,...,p. 



These values are closely related to the so-called 't Hooft parameters in the gauge theories. 



2.2 Critical measures on the real line 

We can make additional progress under the assumption that the external field is analytic on 
the real line. More exactly, we assume that there is a domain Q containing R and a function 
<3? € H(Q) such that cp(x) = <&(x) = Re$(x), x G K. 

Theorem 2.2. With the assumptions and notation introduced above, there exists an analytic 
function R = Rt in real-valued on M, such that 

(c x *(z) + &(z)Y = R(z), zen\S t , (21) 
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where 

c\z)= /^M 

J y - z 

is the Cauchy transform of Xt . 



Formula (21) can be obtained considering variations of the plane of the form x — > x T = 
x + rRe 7^., with 9 G R, r > and z G C \ R, which induce variations on measures, and 
equating to zero the corresponding first variation of the total energy. In this way, the theorem 
is a direct corollary of a more general theorem on equilibrium measures of S-curves from |4CH 
p. 333] or [58j, and the fact that the real line is an 5-curve in any real-analytic field. The 
notion of the S-property turned out to be crucial in the extension of the results on the real 
line to the equilibrium in a harmonic external field in C, considering the complex plane as a 



conductor. In the described situation the representation (21) remains valid, as well as many 
of its corollaries derived in this paper, with some obvious modifications. However, we will 
not go into details here and refer the reader to [58] for a more thorough discussion. 



The proof of Theorem 2.2 actually carries over to a more general situation, e.g. when <j?' 
is meromorphic (in this case R might have poles at the location of poles of <£'). In a more 
physical context, ( |21[ ) was obtained for instance as the saddle-point approximation to the 
large N limit of the number of planar Feynman diagrams of quantum field theory [20J . The 
meromorphic $' has found applications also in the study of the asymptotic distribution of 
zeros of Heine-Stieltjes polynomials |58J or in the analysis of the phase transition in matrix 
model with logarithmic action |50j . 



Let us point out briefly some corollaries of Theorem 2.2 



Clearly, (21) allows us to recover easily the density of the equilibrium measure Xt G Mt 
on R: 

dX t {y) = -^\R t (y)\dy, y G supp(A t ) , (22) 

and its total potential, 

W^(z) = Re f Z s/R^dy (23) 

with an appropriate choice of the lower limit of integration and of the branch of the square 
root. It follows further (as it was shown in [26]) that the support St of consists of a finite 
number of intervals, whose endpoints are simple real zeros of Rt] the unknown zeros of Rt 
become the main parameters determining Xt and its support. 



In this representation we do not have equations (17) anymore (they are now "embedded" 



in the representation (22), which automatically defines a bounded density). Equations (20) 
now take the form 

Q 2fc+1 

y/R t (y)dy =0, fc=l,...,p-l. (24) 

a 2k 

We could continue rewriting further the conditions on the unknown parameters, taking into 



account e.g. the explicit form of and identity (21 ), but in no case the equivalent equations 



will be elementary, or at least suitable for a truly constructive solution. 
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One approach that renders an effective procedure for investigation of those equations is 
to consider the dependence of the system from its parameters. The main parameter is t, and 
it is the central object in this paper. Details are presented in subsections 2.3 and 3.2 below. 



2.3 A dynamical system 

A "dynamical" description of the support supp(Af) was proposed by Buyarov and Rakhma- 
nov in [21] (see also |7H formula (1.8) and Lemma 1], where the particular one-cut case was 
analyzed). For an external field tp, t > and a Borel measure a G Mt let us define two sets, 

S t (a) : = supp(a) J S\a) := {x G M : W%{x) = imnW°(s)}; 

when a = Xt, we omit the reference to the measure, writing St instead of St(Xt) and *S* instead 
of S t (Xt). Then the equilibrium condition Q can be expressed in the following equivalent 
form: 

S t C S\ (25) 
and Xt is the unique measure in M. t satisfying (25). In the regular case, St = S*. Otherwise, 



points from S t \ St are called in \27\ I5T] singular points of type I, and represent also a "birth 
of a cut" : 

s* = p| S T 

T>t 

(see \21[ Theorem 2]); moreover, for an analytic and real- valued ip, S l \ St is at most denu- 
merable (thus, cap(5** \ St) = 0). The following second main ingredient of our analysis is 



actually a simplified version of the result in [21], but can be derived also from Theorem 2.2 



Theorem 2.3. For the given real- analytic external field ip on M., the family of the equilibrium 
measures Xt is monotonically increasing, and for every t > 0, 

dX t dc t , , . . dl^Xt) , . 

■^=«a. = I[u Si ] = p(S t ), ~^ = 2ct. (26) 

Here lok and p(K) = — logcap(X) are correspondingly the Robin measure and the Robin 
constant of the compact K, and ct is the extremal constant defined in Q . 



In the terminology of the random matrix ensembles, Theorem 2.3 states that, up to a 
factor 2, the first variation of the infinite volume free energy is the equilibrium constant, 
while the second variation is the Robin constant of the support (limit spectrum) St- 



Formulas (26) have been rediscovered in less generality in several works, see e.g. [91 134] . 



The first formula in (26) should be used in general in a weak sense: for any continuous 
and compactly supported function / on R, 

d_ 
dt 



f(x)dX t (x) = / f[x)(ko St (x) 



However, since in our context the equilibrium measures are absolutely continuous with respect 
to the Lebesgue measure on IR, the differentiation of the measure Xt with respect to t in 



(26) can be understood in the strong sense of the standard differentiation of its density 



X't(x) = dXt(x)/dx. 
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Proof. The detailed proof of Theorem |2 . 3| (under much less restrictive assumptions on ip) was 
given in [21]. Here we just point out that for the real-analytic external field formulas (26) 
can be established directly from Theorem 2.2 Indeed, a variation with respect to t of the 
equilibrium property ([8]) and of the normalization condition At(R) = t yield that dXt/dt is a 
probability measure with the logarithmic potential constant on St- Formula (22) shows that 
it is absolutely continuous with respect to the Lebesgue measure on M, and that its density 
is of the form 



1 

2^ 



P(x) 



'Rt{x) 



where P is a polynomial. These property characterize the Robin measure of St, which estab- 
lishes the first formula in (26). Moreover, the constant value of y dXt / dt on St must be the 
Robin constant of St, so that by differentiating ^ we get that dct/dt = p(St). 
Finally, by (§, 



=tc t + J ipdX t . 

From the first two identities in (1261) , just established, it follows that 



dl v [\t 
dT~ 



c t + t- 



.dct 
dt 



d 
dt 



J ipdX t = c t + tp(St) + J 



tpdust, 



and it remains to observe that 



tp(St) + J<pdu St = J V" 3 * dX t + f tpdu St = f {v Xt + if) du St = c t . 



□ 



Remark 2.1. Since ct(<p + s) = ct(y>) + s for s 6 I, we can always assume that the external 
field <p is normalized by 

min (p(x) = 0. 

In \21\ Theorem 2] it was additionally proved that with this normalization, the external field 
ip has the form 

poo 

f(x)= gs T (x,oo)dT, xeMSt, (27) 
Jo t>o 

where gs T (x, oo) is the Green function of the domain C\St with pole at infinity, and the 
extremal constant q has a limit when t — > 0: 



liniQ = mmip(x) = 0. 

t->0 x<=K 



2.4 Variations of the external field 

Here we bring to consideration another useful technique related to the variation of the external 
field. Let <p(x, r) be a one-parametric family of external fields depending on a parameter r. 
Without loss of generality, let r S {—£, e); for simplicity, we assume here that ip is real-analytic 
in both variables. The following result, extending the previous theorem, holds: 
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Theorem 2.4. Let t > be fixed. Assume that (p(x,r) : K x (— e, e) — )■ M is real-analytic in 
both variables, and 

d(p(x,r) 



dr 



(28) 



r=0 



uniformly in x in a neighborhood of St = supp(Ai((/?(-, 0))). 
T/ien 



Mi 



r=0 



where the (in general) signed measure /i is uniquely determined by the following conditions: 

dc t 



supp( M ) C S t , m(^) = 0, + ip 



dr 



const on S^. 



T=0 



Proof. Let /Lt be any limit point (in the weak topology) of the family 

A t (r) - A t (0) 



(29) 



(30) 



(which exists due to its weak compactness). The first property of jx in (29) is obvious; the 
second one is a consequence that measures in ( 30 ) are neutral (their mass is 0) . Finally, with 
our analyticity assumptions on (p it easily follows from (21) (see also |58} Lemma 5.3]) that 

C Xt( - T \x) + ^(x,T) = 0, 

in the interior of the support of At(r). Consequently, 

C»(x)+i>'(x) = 0, xeS t , 

which implies the third equality in (1291) . 



Conditions (29) determine \i uniquely. Indeed, the logarithmic potential of the difference 
of two signed measures [i satisfying (29) is constant and equal to on St. It remains to use 
that the logarithmic energy is a positive definite functional on neutral signed measures (see 
e.g [521 Theorem 1.16]). □ 



Remark 2.2. Clearly, Theorem 2.4 can be extended to more general (e.g., smooth) external 
fields ip, but we prefer to keep the settings as simple as possible, and still sufficient for our 
purposes. 



Remark 2.3. As it was mentioned above (see (12)), we can reduce our analysis to unit (prob- 
ability) measures on M, regarding the total mass or temperature t as a particular parameter 
of the external fields of the form ip/t. Consequently, it is convenient to clarify the relation 
between Theorems 12.31 and 12. 41 



According to (12), 



Xt(ip) = t\l, where A* = Ai (</?(•, t)), <p(-,t) = - ip. 
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The derivative of </?(•, i) at a given t > satisfies (28) with ip = —ip/t 2 , and from Theorem 



2.4 



it follows that dX*/dt is a signed measure characterized by the equilibrium conditions 



(29). In this particular case p can be computed explicitly: observe that 



at 



so that applying Theorem 2.3 we get 



~dt 



X\+t 



1 



dX\ 
dt ' 



(31) 



M = t (ws t - A^) ; 



obviously, the right hand side is a neutral measure (f-t(St) 
equilibrium constant equal to 

d fct\ p(S t ) _ ct 
t 2 ' 



0) satisfying (29), with the 



dt \ t 



t 



The reciprocal statement is also valid: formula (31) shows that the logarithmic potential 



of dXt/dt is constant on St, fact that combined with the monotonicity of Xt yields the first 



formula in (26). 



3 A polynomial external field 

After the brief review of the necessary technical background we have carried out in the 
previous section, we can start our study of equilibrium measures in a polynomial external 
field. We define 



2 m 



<p(x) = tj x j , tj £ R, t 2m > 0, 



(32) 



where the coefficients tj (the coupling constants) are the parameters of the problem. For 
convenience, we adopt the normalization 



t-2i, 



1 

2m ' 



(33) 



as well as include the total mass or temperature to = t as an additional coupling constant. 
In this way, the equilibrium measure A = \t(ip,M) is a function of 2m real parameters tj, 
i = 0,l,...,2m-l. 

3.1 /^-representation and (A, £> ^representation 

First, we reformulate Theorem |2.2| and its corollaries specifically for the case at hand, that 



is, when <p is defined by (32)-(33). The theorem implies that for the associated equilibrium 



measure At G Mt with total mass t > and support St C M, the function 



R t (z) = [C Xt + <p') (z) = z im ~ 2 + lower degree terms 



(34) 
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is a polynomial in z of degree Am— 2 with coefficients depending of t. Actually, if we denote by 
(/)© the polynomial part of the expansion of / at infinity, then the normalization A^(M) = t 
implies that 



Rt = {iff + 2(<p'C 



[if') 2 (z) — 2t z 2m 2 + lower degree terms. 



(35) 



Let a\ < CL2 < ■ ■ ■ < ci2p, 1 < p < m, denote the real zeros of Rt of odd multiplicities. 
Observe that Rt cannot have complex zeros of odd multiplicity, and each real zero of odd 
multiplicity of Rt is an end point of St, so that expression (15) holds. We define monic 
polynomials A and B by 



A(z) 



2p 
k=l 



dk), and B 2 (z) 



Rt(z) 
A{z) 



(notice that Rt(z)/A(z) is a monic polynomial of degree 2{2m — p — 1) whose zeros have even 
multiplicities and, therefore, the monic polynomial B of degree 2m — p — 1 > is correctly 



defined). By (24), for p > 1 there should be at least one zero of B in each gap [a2k, a2fc+i]; 
this readily yields the bound p < m. Naturally, both polynomials A and B, as well as the 
endpoints at of the support St, are functions of the parameter t, fact that we usually omit 
from notation for the sake of brevity. 

In what follows we understand by A 1 / 2 and R\^ 2 the holomorphic branches of these 
functions in C \ St determined by the condition A l / 2 {x) > and R 1 J 2 {x) > for x > a2 P , as 
well as /+ denotes the boundary value of / on K from the upper half-plane. 



By (22) 



d\ t (x) 



Til 



(i?y 2 ) + (x) dx, x G St- 



(36) 



The factorization 



Rt(z) = A{z)B 2 (z) = A(z;t)B 2 (z;t) 



(37) 

plays a fundamental role in the sequel. In particular, the ^-representation of the the 



equilibrium measure in (36) may be equivalently written in terms of the following (A,B)- 
representation: 

d\ t ( x ) = — A 1 / 2 (x)B{x) dx, x e S t . (38) 



TTI 



In a one-cut case (p = 1) the equalities (35) and (37) render a system of 2m algebraic 
equations on the zeros of A and B, matching the number of unknowns. Taking into account 
additionally that for p = 1 equation (17) is algebraic in a^-'s, this shows that in this (and ba- 
sically, only in this) case all parameters of the equilibrium measure Aj are algebraic functions 



of the coefficients of the external field (32) 



We conclude this section with two trivial but useful identities: 



A^B 



R 1 ' 2 



B 



A 1 / 2 



and 



(39) 



(40) 
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(see (35}). 

The (A, ^-representation introduced above has several advantages: it allows us to recast 
the discussion on the evolution of the support St and of the phase transitions for the free 
energy in terms of the zeros of A and B. This will be carried out in the next two subsections. 
A more general analysis of the variation of St as a function of all coupling constants is done 
in Subsection 13.41 



3.2 A dynamical system revisited 



Our first task is to rewrite the differentiation formulas of Theorem 12.31 in terms of the zeros 
of A and B, which can be regarded as phase coordinates of a dynamical system. 

Recall that if St has the form (15), then its Robin measure uj = oos t (see the definition in 

h(x) 



Section 1.2) is 



1 h(x) 1 

dust [x) = — —rp, dx = - 

m A^ 2 {x) * 



dx, x £ St, 



while for the complex Green function G(z) = Gs t (z; oo) with pole at infinity we have 

h(z) 



G'(z) 



zeC\S t . 



(41) 



Here h is a real monic polynomial of degree p — 1; for p = 1, h = 1, while for p > 1 it is 
uniquely determined by conditions 



a 2k + l 



h[x) 



a 2k 



dx = 0, k = 1, . . . ,p — 1. 



It follows in particular that for p > 2, 



p-i 



K x ) = Y[( x ~ Cj), Cj e [a 2j ,a2j+i], j = l,...,p-l. 

i=i 

In order to denote the differentiation with respect to the total mass t we will use either 
d/dt or the dot over the function, indistinctly. Recall that at are the zeros of A, and denote 
by bk the zeros of B. 

Theorem 3.1. If A and B have no common zeros then 

2h(a k ) 



a k 



b k = Res 



A'(a k )B(a k ) 
h(x) 



x=b k A(x)B(x) 



k = l,...,2p, 
, k = 1, . . . , deg(B) = 2m — p — 1. 



(42) 



If additionally all zeros of B are simple, then the second set of equations in (42) simplifies to 
; h(b k ) 



A(b k )B'(b k ) 



1, . . . , deg(-B) = 2m — p — 1. 
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Proof. Taking the square root in both sides of ( 34 ) and differentiating the resulting formula 
with respect to t (recall that ip does not depend on t), we get 



Rt{z) 
2R\'\z 



d 
dt 



By Theorem 2.3 and (41), 
d 



dt 



C Xt (z) = C"{z) = -G'(z) 



C Xt (z). 



Hz) 



and further, with the help of (37), 



BA + 2AB 



-2h 



or 



AV2(*)' 



A B 
A +2 B 



zeC\St, 



2h 
AB' 



(43) 



Assertions of the theorem follow by equating residues of the rational functions in the last 
equation above. □ 



Remark 3.1. A subset of equations from (42) corresponding to dk has been obtained in several 



places before, in particular, in the relevant work of Bleher and Eynard [9], although using a 
different approach. That paper is, also probably, one of the first works where some analytic 
properties of phase transitions in Dyson gases were rigorously studied. It contains several 
noteworthy results, such as a discussion of the string equations and the double scaling limit of 
the correlation functions when simultaneously the volume goes to infinity and the parameter 
t approaches its critical value with an appropriate speed. 

Equations (42) are similar in form to a system of ODEs studied by Dubrovin in [30J 



for the dynamics of the Korteweg-de Vries equation in the class of finite-zone or finite-band 
potentials. Curiously, Dubrovin's equations govern the evolution of the "spurious" poles 
of the diagonal Pade approximants to rational modifications of a Markov function (Cauchy 
transform C M ) [81) . and are equivalent to the equations obtained by one of the authors in 
|68j in terms of the harmonic measure of the support of ji. 



3.3 Classification of singularities 



We will see in the sequel that conditions of Theorem 3.1 are satisfied for any t except for a 
finite number of values. These values of t may be called critical; they correspond to some of 
the phase transitions but not to all of them. There are other significant values, also presenting 
phase transitions, which are not critical in the above mentioned sense. 

As it follows from Theorems 2.3 and |2.4| (see also [13] and |5H Theorem 1.3. (iii)]), the 
endpoints of the support of the equilibrium measure, its density function, and the corre- 
sponding equilibrium energy (or the infinite volume free energy) are analytic functions of the 
coefficients of the external field, except for a finite number of values where the analyticity 
breaks down. These critical points divide the parameter space into different phases of the 
model and represent phase transitions. Some of them (but not all) correspond also to the 
change of topology of the support of the equilibrium measure. 
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Recall the classification of the singularities of the equilibrium measure we summarized in 



Subsection 1.4 From the point of view of the evolution in t and with the R- (or (A,B)-) 



representation at hand these singularities are now generically classified as follows. 

• Singularity of type I is a bifurcation b — > (a + ,a~), representing a birth of the cut, 
is the event at a critical time t = T when a simple zero of B ceases to exist and at its 
place two new zeros of A are born. Formally, at t = T the inequality Q, 

W^{x) = V Xt {x) + <p(x) > max^(i), 

^ x&St 

is no longer strict in E \ St, and the equality is attained at some point 6 where 
we will have B(b) = 0. At this point a bifurcation of the zero b occurs. 

A significant property of this situation is that the phase transition occurs by satura- 
tion of the inequality ([8]); the moment t = T of the bifurcation is not defined by the 
dynamical system, i.e. it is not its singular point. All the phase parameters cij,bk may 
be analytically continued through t = T. The solution of the system for t > T would 
give us a critical, but not an equilibrium, measure. 

• Singularity of type II is the opposite event (a + ,a~) — > b, or fusion of two cuts, 
consisting in the collision and subsequent disappearance of two zeros of A (note that a 
zero of the complex Green function, "trapped" in the closing gap, disappears simulta- 
neously), and an appearance of a (double) zero of B, followed by the splitting of this 
real double zero into two complex simple zeros. The collision of two zeros of A is a 
critical point of the dynamical system. 

We note that as a rare event (event of a higher co-dimension) it may happen that a 
number of other cuts were present in the vanishing gap (a~,a + ) immediately before 
the collision; they all disappear at the moment of collision of a~ and a + . This will be 
accompanied by an appearance for a moment of a zero of B of an even multiplicity 
higher than two. 

Thus, a double zero (or of an even higher multiplicity) of B is present at the moment 



of collision inside St- According to (38), the density of At vanishes at these points with 
an integer even order. 

Singularity of type III are the endpoints of St where Rt has a multiple zero; they 



correspond to the case when A and B in (37) have a common real zero. 



Additionally, a special situation is created when two complex-conjugate simple zeros of B 
collide on M\St, and either bounce back to the complex plane or continue their evolution as 
two real simple roots of B (at this moment, two new local extrema of the total potential W Xt 
in R \ St are born). At these values of the parameters the free energy is still analytic, so this 
is not a phase transition in the sense we agreed to use in this study, although the colliding 
zeros of B lose analyticity with respect to t. In a certain sense, the singularity of type III is 
a limiting case of this phenomenon, when it occurs at an endpoint of St- 

The reader should be aware of a certain freedom in our terminology regarding the sin- 
gular points: we refer to singularities meaning both the value of the parameter t at which 
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a bifurcation occurs, and the point on M where the actual bifurcation takes place. We hope 
that the correct meaning in each case is clear from the context and will not lead the reader 
to confusion. 

We return to the analysis of the local behavior at the singularities of the system in terms 
of t = to in Section |4| 



3.4 Other variations of the external field 

So far we have been regarding the equilibrium measure At and its R- and (A, ^-representations 
as functions of the total mass t, assuming that (p is fixed. Now we discuss a more general 
problem: the dependence of A and, correspondingly, of R, A, and B from all the coefficients 



tj of the external field (32). It is a remarkable fact that the differentiation formulas with 
respect to the coupling constants tj, j = 1, 2, ... , have the same form as the differentiation 
formula with respect to t = to, and that they all can be obtained in a unified way. 
We start with a simple technical observation: 

Lemma 3.1. Given points a\ < ai < ■ ■ ■ < a<i v and the corresponding monic polynomial 
A{z) = n)c=i i z ~ a k), there exist polynomials hj, deg(hj) = j + p — 1, j = 0, 1, . . . ; uniquely 
determined by the following conditions: 

ho (z) = -zP- 1 + ..., ^ z) =^ J ~ 1 + °(^)^ z ^°°> 3>h (44) 



and 

AVHx 



z- 



a 2fc+i h(x) 

n ' dx = 0, k = l,...,p— 1, j = 0,l,.... (45) 



'1 2k 



Proof. Clearly, ho(x) = —h(x), where h is the numerator of G' introduced in Subsection 3.2 



see (41). 



Furthermore, for each j £ N, (44) means that the coefficients corresponding to powers 
z~ l ,1, z, . . . , zi~ 2 of the Laurent expansion of hjA~ x l 2 at infinity vanish, which renders j 
linear equations, additional to p — 1 linear equations (45) on the p + j — 1 coefficients of hj. 



The corresponding homogeneous linear equations are obtained by setting 



A 1 / 2 {z) 



-=• , z — y oo, 



along with (45). In particular, every hj is of degree at most p — 2, and according to (45), has 
a zero in each interval (a2fc, a,2k+l), k = 1, . . . ,p — 1, which yields only the trivial solution for 
this system. □ 

Since we are going to write all differentiation formulas in a unified way, we prefer to use 
here the notation 

d 

d i '■= QjT> j = 0, 1, ... ,2m - 1, 

with to = t. 
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Theorem 3.2. Let the polynomial external field p be given by (32) -(33). Let also Xt = 
\t(ip,W) £ Ait denote the corresponding equilibrium measure of mass t > 0, and let polyno- 

(46) 



mials A and B be the (A, B) -representation of this equilibrium measure, see (38). 
Then 



Ay\z\ 

where polynomials hj are given in Lemma \3A 
Moreover, 

2m- 1 



zeC\S t , j = 0,l,...,2m-l, 



tj dj (C Xt + if') (z) = [C Xt + v?') (z), zeC\S t . 

j=0 



In consequence, 



and 



B djA + 2AdjB = 2hj, j = 0, 1, . . . , 2m - 1, 

2m- 1 

^ t j h j (z) = A(z)B(z). 

3=0 



Proof. For j = 0, formulas (46) and (48) are just a restatement of Theorem 3.1 



Furthermore, from Theorem |2.4| it follows that for j = 1, . . . , 2m — 1, 

djC Xt (z) = C 9 J Xt (z) = C^(z), 
where uij is a signed measure on St satisfying 

LOj(St) = 0, (x) + x j = dj = const on S t . 



It follows from (44)-(45) that the multivalued analytic function 

hj(x) 



Gj{z) 



■ dx, z G C \ St 



'■MP A1/2 ( X ) 

has a single- valued real part, which is continuous in C and satisfies 
Re Gj{z) = 0, x G S t , and Gj(z) = z j + O (1) , 



Z —f oo. 



Taking into account (50) we conclude that 



ReGj(z) = V"i(z) + z j 



z e 



Since for j = 1, . . . , 2m — 1, 



djf(x) = x 3 , djif'(x) = jx- 



(47) 

(48) 
(49) 



(50) 



identity of (46) follows for all remaining j's. 

Finally, Wm shows that C A * + <p' is a homogeneous function of degree 1 of the vector of 
coupling constants to, t\, . . . , so that (47) is just Euler's theorem for such a function. Formula 
(49) is obtained by replacing (34) and (46) correspondingly in the left and right hand sides 
of (47). □ 
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Evaluating (49) at the zeros of Rt we obtain a set of algebraic identities, 



2m 



^2tj hj(a k ) = 0, k = l,...,2p, 
3=0 
2m 



(51) 



3=0 



0. 



B(z)=0 



called hodograph equations in [lj. Solving them we could find the main parameters of the 
equilibrium measure and of its support. However, their explicit character is misleading: in 
the multi-cut case the dependence of the coefficients of hj from the coupling constants tj's is 
highly transcendental, and as the authors of pQ point out, except for the simplest examples, 



equations (51) are extremely difficult to solve, "even by numerical methods". 



Remark 3.2. Alternatively, following the general methodology put forward in |49| . we can 
derive the identities on the endpoints of the connected components of St considering the 
hyperelliptic Riemann surface 1Z of w 2 = A(z) and its deformations depending on the set 
of coupling constants tj, imposing the condition that the partial derivatives with respect to 
the parameters tj of the corresponding normalized Abelian differentials of this surface are 
given by a meromorphic differential T> on 1Z. This is equivalent to the set of the so-called 
Whitham equations [H] on 1Z. Actually, polynomials hj defined in Lemma 3.1, appear in 
the explicit representation of these normalized Abelian differentials of the third (j = 0) and 
second kind (J £ N) on 1Z. The key connection with the equilibrium problem is provided by 
identity (34), which shows that the differential (ip (z) + R t (z))dz, with tp given by (32), 
can be extended as this meromorphic differential T> on 1Z. This approach was used in [22] to 
obtain in particular an analogue of ( 46 ) , and developed further in pQ . 



Again, a direct consequence of Theorem |3.2| is the possibility to rewrite the differentiation 



formulas (48) in terms of the zeros at of A and b k of B: 



Theorem 3.3. If under assumptions of Theorem 3.2, A and B have no common zeros then 
for j = 0,1,. . .,2m- I, 



2h 3 {a k ) 

U 'j a k — 777 \ T-y / 7j ft — l,...,Zp, 

A'(a k )B(a k ) 



hj(x) 



(52) 



x=b k A(x)B(x) 



, fc = l } ...,deg(B). 



If additionally all zeros of B are simple, then the second set of equations in (52) simplifies to 

hj(bk) 



djh 



A(b k )B'{b k ) 



, fe = l,...,deg(B). 



Proof. It is a consequence of ( 48 ) that 



diA 



(x) + 2 



B 



[x 



AB 



(x), i = 0,l,...,2m-l. 



(53) 
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Observe that 



d A 2p 2p d a 

i -(x) = djlog(A)(x) = ^2djlog(x - a k ) = 



A ^— ' " ^— ' x — ah 

k=l k=l K 



Analogous formula is valid for djB/B. Hence, both the left and the right hand sides in (53) 
are rational functions in x, with possible poles only at the zeros of A and B. The necessary 
identities are established by comparing the corresponding residues at each pole. 

For instance, with the assumption that the zeros of A and B are disjoint, the residue of 
the left hand side of (53) at x = a k is equal to —dja k , which yields the first set of equations 



in ( |52[ ) (recall that by construction, all zeros of A are simple). The analysis of the residues 
at x = b k gives us the remaining identities. □ 

The proof shows how the statement can be modified in the case of coincidence of some 
zeros of A and B (in other words, in the case of roots of Rt of degree higher than 2). Moreover, 



the evolution of Rt is such that as long as the right-hand sides in (52) remain bounded, all 
zeros of and hence, Rt itself, are in C . 



We can rewrite the equations on a k s in (52) in a weaker form: 

hi(a k )dja k = hj(a k )dia k , k = l,...,2p, i, j = 0, 1, . . . , 2m - 1, 
which are the Whitham equations in hydrodynamic form (see [TJ Eq. (77)]). 

Example 3.1. The simplest case to consider is when A(x) = (x — a\)(x — 02), so that St 
consists of a single interval [01,02]. In this situation, ho(x) = — 1, h\(x) = x — (ai + ct2)/2, and 



under assumptions that A and B have disjoint zeros we obtain from (52) that for j = 1,2 



d aj = ——. — ^— — r and <9ia 7 - = — „ . . 
3 A'{a j )B{a j ) 3 B(a,j) 



This can be rewritten as a system of PDE, 

Oiai d ai = 0, 

a .02-01 

Oia 2 H o a 2 = 0, 

which is a rescaled form of the continuum limit of the Toda lattice (see [5]). This set of 
equations for the endpoints of St in the one-cut case was derived previously in [3 H] . 
In an analogous fashion, 

. . 2 / \ (01 — 02)^ 
h2(x) = 2x -(a 1 + a2)x , 

3 3 3 

h 3 (x) = 3x 3 - -(ai + a 2 )x 2 - -(d - a 2 ) 2 x - — (ai - a 2 ) 2 (ai + a 2 ), 
2 a lb 
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which yields the following differential relations: 

3ai + a 2 „ 3(4a? + (ai +a 2 ) 2 ) 

a 1 + 3a 2 and 3(4a| + (o 1 + a 2 ) 2 ) (54) 

9202 " " zbkT' ^ 2 = • 

Again, these can be rewritten as dynamical systems, 

(a 2 - a 1 )(3a 1 + a 2 ) 
9 2 ai 9 ai = 0, 

„ , (02 - oi)(oi + 3a 2 ) 

(? 2 a 2 H d a 2 = 0, 

and 

3(a 2 -ai)(4a 2 + (ai+a 2 ) 2 ) 

03 a 1 — d ai = 0, 

lb 

83O2+ 3(a 2 - O1 )(4 o ; + (» 1 + O2 ff 6ba2 = 
lb 



(55) 



4 Local behavior at phase transitions in the polynomial ex- 
ternal field 

In this section we take a closer look at the behavior of the equilibrium energy for the external 



field (32)-(33) when these phase transitions occur, considering only the variation of the total 



mass (equivalently, temperature or time) t. In order to simplify our analysis, we make some 



generic assumptions on the three types of phase transitions described in Section 3.3 Along 
this section we denote by t = T the critical time at which the phase transition occurs, and 
use the R- and (A, .^-representations of the equilibrium measure Ai as introduced in Section 
IO 

• Singularity of type I: at a time t = T a real zero b of B (a double zero of Rt) splits 
into two simple zeros a_ < a+, and the interval [a_,a+] becomes part of St (birth of a 
cut); we assume that for t in a neighborhood of t = T, b is a simple zero of B; 

• Singularity of type II: at a time t = T two simple zeros a 2s and a 2s +i of A (simple 
zeros of Rt) collide (fusion of two cuts). 

• Singularity of type III: at a time t = T a pair of complex conjugate zeros b and b of 
B (double zeros of Rt) collide with a simple zero a of A, so that \' T (x) = 0(\x — a| 5 / 2 ) 
as x — > a. 



Additionally, in Subsection 4.4 we analyze the scenario when at a time t = T a pair of 
complex conjugate zeros b and b of B (double zeros of R±) collide at bo £ R \ St and either 
bounce back to the complex plane or become two simple real zeros 6_ < b + of B (these real 
zeros are new local extrema of the total potential W Xt on K \ St). 
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Obviously, for a general polynomial potential some of these phase transitions can occur 
simultaneously: for a given value t = T two or more cuts could merge, while a new cut is 
open elsewhere, together with a type III singularity at some endpoints of St- Still, the basic 
"building blocks" of all phase transitions are precisely the four cases described above, which 
we proceed to study. 

It was established in [13J that in the second case, when at t = T two components of St 
merge into a single cut in such a way that ip is regular for t £ (T — e,T), the equilibrium 
energy can be analytically continued through t = T from both sides. For the case of a quartic 
potential cp it was shown in [9] that the energy and its first two derivatives are continuous 
at t = T, but the third derivative has a finite jump. This is a third order phase transition, 
observed also in a circular ensemble of random matrices [H] . With respect to the singularity 
of type III, it is mentioned in |12| that in this case the free energy is expected to have an 
algebraic singularity at t = T, "but this problem has not been studied yet in details" . 

In this section we extend this result to a general multi-cut case and show that for the three 
types of singularities the energy and its first two derivatives (but in general not the third 
one) are continuous at the critical value t = T. The character of the discontinuities is also 
analyzed, and we can summarize our findings as follows: in all cases there is a parameter, that 
we call 5, expressing geometrically the "distance" to the singularity. In the case of a birth of a 
new cut, this is the size of this new component of St; for a fusion of two cuts, this is the size of 
the vanishing gap, while for the singularity of type III we can take it as the distance between 
two colliding zeros of Rt- In the three cases the first two derivatives of the equilibrium energy 
/^[Af] are continuous, while the third derivative has a discontinuity Ap(St). We gather in the 
table below a rough information about the character of the dependence of these parameters 
from At = \t — T\ at the corresponding singularities, see Subsections 4.1-4.3 for a detailed 
asymptotics and further discussion. 





Type I 


Type II 


Type III 


6 


0(y/At/]ag(At)) 


0(VAi) 


0((At)V3) 


Ap(5 t ) 


0((Atlog 2 (At))- 1 ) 


0(1) 


©((At)" 2 / 3 ) 



Clearly, all zeros a^, bk of Rt (and hence, Rt itself) are continuous functions oft. Moreover, 



it follows from (42) that when a singularity occurs, all a^'s and b^s not involved in the phase 
transition are smooth functions of t, with the degree of smoothness, roughly speaking, 1 more 
than the smoothness of 5 in the corresponding column of the table above. 



4.1 Singularity of type I: birth of a cut 

Recall that we assume that at t = T a simple real zero b of B (a double zero of Rt) splits 
into two simple zeros a_ < a + , and the interval [a_, a+] becomes part of St- 

First of all we want to estimate the size of the new cut as a function of t. We use the 
notation introduced above, indicating explicitly the dependence from the parameter t. Let 
us remind the reader in particular that polynomial h is the numerator of the derivative of 



the complex Green function, see (41 ). Our assumptions can be written as follows: for a small 
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e > 0, there exist polynomials A, B and h, such that 
A(x;t) -- 



B(x;t) 
h(x; t) 



A(x;t), te(T-e,T], 
(x - a-)(x - a + )A{x;t), te(T,T + e), 

(x -b)B(x;t), te(T- e,T], 
B{x;t), t€(T,T + e), 

h(x;t), te(T-e,T], 
(x-()h(x;t), t€(T,T + e), 



where a±, b and £ are real- valued continuous functions of t such that 

a_(i = T+) = a+(t = T+) = b(t = T-) = ((t = T+); 

we denote this common value by bo; this is the place where the new cut is born at t = T. 
Remember that it is imposed only by the saturation of the inequality constraint in ^ outside 
of St- Polynomials A, B and h are continuous with respect to the parameter t G (T—e, T+e), 
but represent, generally speaking, different real-analytic functions of t for t < T and t > T. 

Let us denote Q(x) = A(x)B(x), omitting from the notation when possible the explicit 
dependence on t. From (42) it follows that 

; h(b) h(b) 

(a+ — a-)Q(a±) {a + — a-)Q{a±) 
It is convenient to introduce two new variables for t > T: 

£ = £(t) = (a+ + a_)/2, 6 = 6(t) = (a+-a-)/2. (57) 
Adding both equations in ( |56[ ) we get 

• If h(a + ) h{aJ) 



a + — a_ \Q{a + ) Q(a_ 

1 /(a+-C)M«+) (a_-C)M«-) 



\ rv , n, ' ' i £ (T, T + e), 

a+ — a_ \ Q{a+) Q{ct-) 

which shows in particular that 

h(b 



Km « - ' *<*> 



and we conclude that the function defined piecewise as 



= lim b, 

x =b Q(Pq) t^T- 



b{t), t<T, 
£(*), t>T, 
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is also C 1 in (T-e,T + e). 

Analogously, subtracting equations in (J56|) yields that for t G (T, T + e) 



A fA 2 ^ = , M^-) = (q+ - C)M«+) , (a- - C)h(a-) 

dt [ ' Q(o+) Q(o_) ' Q(a+) Q(a_) 



or equivalently, 



(it 



= (a_ - C) 



, MO \ . fe(o+) 
Q(a+) Q(0/ Q(a+)' 



(58) 



Let us denote by a the largest zero of A satisfying a < bo. Then condition (45) reads as 

/ \ -i 



x — £ 



dx = 



(c-a-r 1 



dx 



x — a_ 



(59) 



Clearly, 



lim 



x — o_ 



^T+y a yz(xTt) A 7a(xTT) 



6o 



dx 



,dx / 0. 



On the other hand, 

dx 



r 

J a 



dx 



\/A(x;t) y/A(a-;t) J a i/(a_ - x)(a+ - x)) 



+ 



F 

J a 



yftal -x)(a+-x)) y^Afcjt) ^(a-;*)^ 
1 log i + 0(1), <5^0 + . 



dx 



From ( 58 ) and ( 59 ) it follows now that if we denote 

h(b ) f bo dx 



Q(bo)J a y/A(x;T) 



dx > 0, 



then 



d 



(5 2 



-C + oil) +25 h^_ + 



dt"~ ' log(,5 2 ) + 0(l) Q(b , 
Notice that the solution of the ODE f'(x) = a/(log(/(x) + b), with /(0+) = 0, satisfies 

' log log x 



(60) 



fix) 



ax 



log(/(x)) + b - 1 log(x) 



logx 



, x ->• 0+, 



so that (60) implies that 



5 2 (t) = -C 



t - T 
log(t - T) 



l + O 



logiog(t-r) 

log(t - T) 



t^T + 



(61) 
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This is consistent with the log(n) jn scaling used e.g. in 

We turn next to the asymptotic behavior of the Robin constant p(St) of the support St, 
which according to Theorem 2.3 is the second derivative of the infinite volume free energy. 
Let us study first the following model situation, from which the general conclusion is readily 
derived. 

Assume that E C M is a union of a finite number of disjoint real intervals, and A5 = 
[£ ~~ <5) £ + <5] is disjoint with E for 5 < e. We denote Eg = EuAs, and slightly abuse notation, 
writing p{5) instead of p(Es). Observe that with this notation, p(0) = p(E). Our goal is to 
study the behavior of p(5) as 5 — > 0+. 

As usual, gE(z,a) stands for the Green's function of C \ E with a pole at a ^ E. The 
well-known relation holds: 

g E (z,0= log 1^1+7(0 + 0(1), z^t, (62) 

which defines the constant y(Q- 

Lemma 4.1. There exists a constant C > (uniform in E and such that for 5 < e, 

si(f.°°) 



" (4) -" (0) + 7 (0 + M2/ 4 ) 



< C5. 



Proof. Let p be the balayage of the unit point mass at £ onto E, so that (see [52] ) 

^(S) = fog j-i— + OO). (63) 



Comparing it with ( 62 ) we conclude that 

g E (.z,0=g E (Z,oo)-V' t (z)-log\z-t\. (64) 

Indeed, the function in the right hand side is harmonic in C\ (ED {£}), vanishing on E, and 
has the appropriate logarithmic behavior when z — > £. 

Recall that we denote by oj e and be the Robin measures of E and A§, respectively. 
For a parameter m > define 

a = as = oje + muj^ s — mp; (65) 

observe that since p! is bounded on E, and u' E is bounded away from zero on E, a is a 
positive measure for sufficiently small values of m > 0. Consider the potential y CT for such 
values of m; we want to find its bounds on E$. 
Let us start with x £ E. We have 

V a (x) - p(0) = (y" E (x) - p(0)) - m \V^{x) - lo, ' 



+ m [ V^^s (x) - log -. Al 

x -€\. 
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The first parenthesis in the right hand side is = by definition of uj e ; the second one is equal 
to <?.e(£, oo), where we use (63). So, we need to estimate the third parenthesis. 
We have 

V^s (x) = log 



so that 



V WA s (x) - log 



x-C+^ix-O 2 -^ 
2 



log- 



1 + 1/1 



, \x-(\>6. 



Using the straightforward bounds 

2 



log 



1 + Vl - s 2 



log 1 + 



(i + ^/r^2) 2 



we conclude that 



< V"*s (x) - log 



< 



x — £| \x — £ 



\x - £\ > 5. 



Gathering all estimates in (66) we see that 



mgE(C, oo) < V a (x) - p(0) < -mg E (C, oo) + m 



x G E. (67) 



dist{i,E) 

Now we look for similar bounds on A^. Observe that 

V a (0 - P(P) = (V WB (C) - p(0)) - mV (0 + mV UA s (f ). 
By the well-known property of wg, the first parenthesis is equal to — 5^(^,00); by identity 

^(0=^,0°) -7(0, 



PI), 



while by the equilibrium condition (11), V"^* (£) = log(2/<5). Thus, 

y a (0 - P(0) = 00) - m {g E (C, 00) - 7(0) + m log 

Now we set the value of the parameter m to guarantee that 

V°(0-p(0) = -mg E (Z,oo); 

in other words, we take 



m 



gg(e,oo) 

7(0 + log I 



> 0. 



Observe that m — > 0+ for (5 — > 0+, so that we can choose e > such that cr^ defined in (65) 
is a probability measure for < 5 < e. 
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We are ready to establish the sought bounds on A$. Since V UJA s is constant on Ag, we 
have 

V a (x) - V a {£) = {V WE {x) - V" E {£)) - m (V^x) - . 
Clearly, both terms are harmonic on As, so we conclude that 

\V a (x)-V a (0\<C 1 5, 5^0+, (69) 

where C\ is uniform with respect to small variations of E and £. 

It remains to gather (67), (68) and ( |69| ), and use the well-known property 

min V a (x) < p(S) < maxV a (x), 

x£E s x£E s 

which concludes the proof of the lemma. □ 

In order to use this result in the situation of a birth of a cut it is enough to notice that 
Q is a differentiate function of t at t = T, and so is Taking into account (61) we get 

(q , , q v 9 2 s T (bo,oo) 1 / / k2glog(t-r) \\ 
P (S t )-p(S T ) = ^t^ ) { 1 + °{ log(t _ T) JJ> ^ T +- 

In particular, jo(5^) is continuous but non-differentiable at t = T, since it has an infinite 
jump: 

% (St) = - *t^> 1 ^ fl + f»^iy «^T +1 (70) 



<ir v 2 (t - T) log 2 (t — T)\ V log(t - T) 

while the limit 

lim —p(St) 

obviously exists and is finite. 

Expression (70) is consistent with the O ((t — T) -1 log~ 2 (t — T)) jump found in pQ for 
the quartic potential (the simplest transition from one to two cuts, see Section [5]). 



4.2 Singularity of type II: fusion of two cuts 

We turn our attention to the critical transition when two intervals from St collide and merge 
into a single interval. More precisely, we assume that at a time t = T two simple zeros ct2 S 
and <22 S +i of A (simple zeros of Rt) coalesce. 

Again, we start by estimating the size of the gap between two merging cuts as a function 
oft. 
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Using the notation introduced above, our assumptions can be written as follows: for a 
small e > 0, there exist polynomials A, B and h, such that 



A(x;t) 
B(x;t) 
h(x; t) 



(x - a-)(x - a+)A(x;t), t G (T - s,T), 
A(x;t), t€[T,T + e,T), 

\x - b)B{x;t), te(T-e,T), 
(x - b)(x - b)B(x-t), t£[T,T + s,T), 

\x-C)h(x;t), te(T-e,T], 
h(x;t), t€(T,T + e), 



where a± , b and £ are continuous functions of t such that 

a-(t = T-) = a + {t = T-) = b(t = T±) = = T-); 

again, we denote this common value by bo . Notice that a±, b and £ are real- valued for 
t G (T — e,T], while b G C \ M for t G (T, T + e). Polynomials A, B and /i are continuous 
with respect to the parameter t G (T — e, T + e), but represent, generally speaking, different 
real-analytic functions of t for t < T and t > T. 

We denote Q{x) = A(x)B(x), omitting from the notation when possible the explicit 
dependence on t. From (42) it follows that 

. 2(C-o±) 



a± 



o-)(o± - 6) Q(a±) 



C-6 



(6-o_)(6-a+)0(6) 



1 h(6) 



-1 h(b) 



b-bQ(b)' 



b-bQ(b) 



t G (T,T + e). 



, t & (T — e,T), (71) 



(72) 



Adding the first two equations in (71) and using the notation introduced in (57) we get 

j = -1 / (C-Q+)Ma+) (C-a-)h(a-) \ t e(T-E,T), 
a + — a— \(a+ — b)Q(a + ) (a_ — b)Q(a-) J ' 

so that 



lim £ 



Analogously, from (72), 



(C — x)h(x) Y 
(s - 6)Q(x) 

lim 4 (2 Re 6) 



x=6o 



Q(x)J 

and we conclude that the function defined piecewise as 

few, t<T, 

1 2 Re 6(f), t>T, 



fe(g) 
Q{x) 



x=b 



x=bo 



is C 1 in (T-e,T + e). 
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Furthermore, by (71) 



d (r 

dt 



Condition (45) reads as 



(g+ - C)H a +) (C-q-)M°-) 

(a+-b)Q(a+) (b-a_)Q(a_)' 
^ = dx = 0. 



With the change of variables x = (a+ — a_)y + a_ in the integrand we obtain that 

= p ptM£ dy f A 1 J - 1 



where 



so that 



9i(y,t) 



h((a+ — a_)y + a_) 
A/ A (( a + - a -)2/ + a -)' 



lim = 



uniformly for y £ [0,1]. Thus, 



lim 



C-a_ 



lim 



a+-C 1 



t-±T- a + — a_ t-»T- a + — a_ 



Analogously, by (24), 



/ (x - b)y/(a+ - x)(x - a-)\f A(x;t)B(x;t)dx = 0. 

J a— 



With the same change of variables as above, we get 

b-a- fX ' rl 



a_i_ — a_ 



f vVy(X--y)92(y;t)dy (J^ \/y(i - y)g2(y;t) 



dy 



where now 
so that 
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(y;*) = #((«+ - a -)y + a-)V A (( a + -a_)y + a_ 



lim g 2 (y-t)=B(b )y / A(b ) 

t-¥T- 



uniformly for y £ [0,1]. Consequently, 

b — a_ 

hm 



lim 



a+ — b 1 



t-»T- a_i_ — a_ t->T- a + — a_ 2 



Combining (74) and (75) we see that 

lim ^ 



lim — f = 1. 



t->T- b — a- t->T- a+ — 6 



(73) 



(74) 



(75) 
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Using it in (73) we conclude that 

lim ± (6>) = 2^) < o. 
t^T- dt y ' Q(bo) 

This formula shows that two cuts can come together only at zeros bo of B for which (ABh) (bo) < 
0. Furthermore, we have 

Now we turn to the simplified model problem. 
Assume that p > 2, 

p 
k=i 

where = [a,2k-ii a 2k\ are pairwise-disjoint intervals, and assume that for certain s G 
{l,...,p-l}, 

0-2s = £ - a 2s+l = £ + S, 



Es = (J A A 



where £ is fixed and 5 — > 0+. No other endpoint a k depends on S. By (41), for G(z) = 
Ge s {z; oo), we have 

with 

5) = JJ (z - a k ) = (z - £ - <*)(* - £ + ^)A(z), 
fc=i 

and 

p-i 

= = ( z -(s)h(z;S), 

fc=i 

where all Cfc G [«2fc> a 2fc+i]> for A; = 1, ... ,p — 1; notice that deg(/i) = p — 2. Zeros do 
depend on 5, so we will write Ofc(^) when we want to make it explicit. Observe also that 
Gs(0) = £ and Cfc(0), 6 {1, . . . ,p — 2} \ {s}, are the zeros of the derivative of the complex 
Green's function for Eq- 

Lemma 4.2. With the notation above, 

US) - Cfc(0) = O (8 2 ) , 5^0+, fc = l,...,p-l. (77) 

Moreover, 

C»(<5) - £ = 5 2 + O (5 3 ) , 5^0+, (78) 

with 



K, 



| E^m-J E At)- ( 79 ) 



2 l^Vf-CifO) 2 ^ £-a fc 

\j^s s SJ w fc/2s,2s+l s 
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Proof. Consider first the vector-valued function T : [0,e) — > assigning each value of 

5 to the corresponding vector (Ci(^), • • • , Cp-i(^))> defined by equations (45). It is clearly 
differentiate for any 5 £ (0,e). If we denote 

ra. 2 k+i h(x;5) 

G k (6;(i,...,(p-i) = / , 1/2 / A N dx = ' k = l,...,p-l, 

Ja 2k 1 \ x \ 0) 



then by the implicit function theorem, 



/ dG 1 

1 dCi 



8 ri 



8Gp- 



dCp-i 



dG v 



\ 


( id \ 







(80) 



Consider the s-th row of the matrix in the right-hand side of (|80j) for 5 — > 0+. We have 

r < +5 (x-( s )h(x;6) 



G s (5; Ci, • • • , Cp-i) 



dx. 



U-S y / (x-^-5)(x-C + 5)A(x) 
Since only the numerator of the integrand depends on ^ , it is easy to see that 
dG s 



lim 

5->0+ d( s 



-it — , = = ^ and lim ^^*_ s = for k ^ s. 



5^0+ d(k 



On the other hand, the minor of the matrix in the right-hand side of (80) obtained after 



eliminating the s-th row and column is clearly invertible for 5 — > 0+: it corresponds to the 



system (80) for Eq, and the endpoints of Eq are free ends where no phase transition occurs. 



Hence, expanding the matrix in (80) along its s-th row we conclude that it is invertible 



for small values of < 5 < e, and we can write 



/ §Gi 
1 SCi 



Observe that for j ^ s, 



\ 1 / o r 



05 

while for j = s, 



QQ r a 2j+ i 

' — 



dGp—i dGp—i 

\~9CT~ ••• ~Wp~^/ 



(x - ( s )h(x;5) 



„ 2i ^( x -^-5Y{x-i + 5YA(x) 
Um dG s _ h{H- 0) 



s^o+ 05 ^-A(£) i V^s 1 



0{5), 5^0+, 



ds = 0. 



This proves (77) 



Let us be more precise about the asymptotics of ( s . By (45) 



(x - ( s )h(x;5) 

8 y/5^-( X -0 2 V^Mx~) 



dx = 0. 
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Thus, defining rj = £ s — £ and making in the integral above the appropriate change of variables 
we obtain an expression for r\: 

where 

/M)= ^ + ^) (81) 

Since /(r, 5) = /(0, <5) + /;(0, <5)r + ^ T (0, 5)r 2 /2+er 3 , and defining M 2 = max H < a |/^(t, 
M 3 = max| T |< 5 |/^ t (t,£)|, we get that 

MoS 2 

dT = f(0,S) + e 1 , \ £l \< ~ 





c 


f(r,S) 


7T „ 




VS 2 - r 2 




i: 


rf(r,5) 


7T . 




V5 2 - T 2 



4 ' 

2 16 



so that 



/;(o,5) , 2 



where 



^mf) +es ' (82) 

£2 -^y(M)/(2/(o,^) ) = { ^ 0+ ^ 

/(0,<5) + £i 

We claim that in the asymptotic expression above we can replace /|(0, <5)//(0, <5) by 
/£(0, 0)//(0, 0) + (<5 2 ). This is a direct consequence of (77) and the fact that in /(t,<5), 
only the numerator depends on (5. 

This proves the lemma. □ 

As in the case of the birth of a cut, we study the asymptotics of the Robin constant 
p(6) = p(Eg) as 6 -> 0+. 

Observe that for x ^ £ and ( s - £ = K 6 5 2 + (<5 3 ) , <5 -> 0+, 

' r_C l +(wr^v>-—^P + 0{5 i ), (83) 



with the 0(<5 3 ) term uniform in x. 

Recall that for all 5 > 0, 6) is a monic polynomial of degree p — 2, and by (77), 

<5) - h(z; 0) = (5 2 ) , 5^0+ (84) 

uniformly in the endpoints of Eq. This motivates us to define 

= lim M*; *)-*(*;<>) 

5^0+ 5 2 
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the existence of this limit will be established next. Meanwhile, it is clear that H is a poly- 
nomial of degree at most p — 3 (for p = 2, function f(r,5) in (81) is a constant, so that by 
([82j , H = in this case). 



Assume that p > 3; equations (45) for k ^ s yield: 

x — Q s h[x;8) h(x;0) 

y/(x - | + <5)(x - I - 5) t/Mx) ~ y/A(x) 




Dividing it through by 8 2 , using (78)-(j79|) and considering limit when 8 
following equations on H: 



dx = 0. 

0+ we obtain the 



X 



a 2k 



dx 



a 2k + l 



"2k 



1 



h(x;0) 

Z 2(x-0 2 J7W) 



dx, 



,p-l, k^ s, 



(86) 

where is defined in (79). This renders a system of p — 2 linear equations with p — 2 
unknowns (the coefficients of H), which has a unique solution (as the consideration of the 
corresponding homogeneous system clearly shows), and in particular, establishes the existence 
of the limit in ( 85 ) . 



From (41 ) we have for y > a>2 P , 

ry h(x;5) 



so that 



G(y,E s ) 
G(y,E ) 

p(8)-p(0)-- 



a 2p yJA(x;$) 
y h(x;0) 

a 2p y/A(x;0) 
h(x; 5) 




dx = log(y) + p(S) + o(l), 
da; = log(y) + p(0) + o(l), 

h(x;0) 



+oo, 



+00, 



^Ajx^I) ^A(x^j 

x — ( s h(x;S) h(x;0) 




dx. 



and the integral is convergent for every sufficiently small 5 > 0. Taking into account (77), 
(83) and (84), we conclude that p{5) — p(0) = 0{5 2 ). Hence, dividing the identity above 
through by 5' 2 and using the definition of H we get 



lim 

5->0+ 



5 2 



0-2p 



H(x) 



K f 



1 



e 2(x-0' 



h(x;0) 



dx 



^Mx)' 



(87) 



We summarize this in the following lemma: 
Lemma 4.3. Under the assumptions above, 



p(6)-p(0) = K p 5 2 (l + 0(l)), 5^0+, 



where the constant K p is given by the right hand side of (87), and the polynomial H is 
uniquely defined by the equations (86) for p > 3, or H = for p = 2. 
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that 



Now we go back to the phase transition when two cuts merge. Using formula ( 76 ) we see 

h(bo) 



p(S t ) - p(S [ 



T 



2K n 



-(t-T)(l + (l)), t 



A{b )B(o 0) 

Observe that this expression involves explicitly the point &o where two cuts merged; clearly, 
this is not the case if we consider the limit of dp(St)/dt as t — > T+. We could conclude from 
here that dp(St)/dt at t = T has finite but, in general, different values from the left and from 
the right. In particular, p(St) is not differentiable at t = T. 

Example 4.1. As an illustration, let us consider the case when St has only two cuts that 
merge into a single interval at t = T. Using the notation above, this means that A(x; t) = 
(x — ai)(x — a 2 ), with a\ < a_ < 60 < a + < a 2 for t 6 (T — e,T), h(x;t) = 1, and H = 0. 
Thus, by m, 



8\t) 



A(b )B(b 



-(t-T)(l + o(l)), t^T 



Since St = [01,02] for t 6 (T, T + e), we have jo(St) = log(4/(a2 — a\))- Using (42), 



dt 



p{St 



a 2 — ai 



et 2 — ai (a 2 - ai) 2 \i?(a 2 ) -B(ai) 



+ 



so that 



lim -z-p(St) 
t^T+ dt 



+ 



(02 - ai) 2 V (02 - b ) 2 B(a 2 ; T) (ai - b ) 2 B{ar,T) 



<0, 



where we take the values = aj(t = T). 

Consider now i G (T — e,T); with the notations above, 



1 1 

+ 



£ - ai £ - a 4 



a L > 



1 arccos(#) 



where 



Furthermore, 



ai + ^2 - 2£ 
02 — ai 



e(-l,l). 



dx 



arccos(#) 



1 



, 02 (x-0 2 VMx) (0 2 -l)(a 2 -a x )2 V 
Gathering these identities, setting £ = 60 and taking into account the expression of 6, we get 



lim — p(St) 
t^ T -dt HK ' 



1 



A(6 ;T) 2 B(& ;T) 



< 0. 



From (88) and (89) we can easily compute the finite jump of -n p(St) at t = T. 
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Further simplifications are obtained for the quartic potential, when B = 1. Equations 
above boil down now to: 

5\t) = - 1 (r-t)(l + (l)), t^T- 

{b - ai)(a 2 - b ) 

(this expression matches, after due transformations, formula (2.41) in [9]), and 

,■ d /on 1 v d /on 2 ( 1 1 

hm. —p(S t ) = --, 7-^, lim — p{b t ) = --, 7 rr^ + 



t->T- dt (a x - b ) 2 (a 2 - 6o) 2 ' t~^T+ dt (a 2 - ai) 2 \ ( a 2 - &o) 2 (oi - 6o) : " 

so that 

lim A - hm £ = - f -. a ; + a2 ~^° -J) 2 < 0. (90) 

t^T+dt' y ' t^T-dt Hy } \(a 2 - a^iax - b Q )(a 2 - b ) J ~ 

In order to compare these formulas with those obtained in [9], we must set a\ = —2 = —a 2 
and bo = 2c\ = 2cos7re, and denote s\ = sin7re. It yields 

lim ~ P(St) ~ Hm 4 P(S t ) = — — A < 0. (91) 
t^T+ dt PV ; t-*T- dt PV ; 16s| ~ v ; 

This value differs slightly from the one in [9], probably due to a minor error in formula (2.55) 
therein. 

Observe finally that that the jump in (90) or (91) is strictly negative, unless bo (place of 
collision) coincides with the midpoint of the interval [01,02]. For the quartic external field 
it takes place if and only if it is symmetric, i.e. attains its global minimum at two distinct 
points. This conclusion is a straightforward consequence of the formulas in (39). Indeed, if 
a\ = —a 2 and 60 = 0, then the first formula in (39) gives us that 

If \ 3 a l / \ 1 2/ 2 2\ 

(p(x) = x ~y x ' f{ x ) = ^ x ( x - a i)- 

Since the situation is invariant by translation in M, we can conclude that if bo is the midpoint 
of the interval [01,02] then tp has two global minima, situated at 

«i + a 2 a 2 - a\ 
2 2^2 ' 

and a global maximum at (ai + a 2 )/2. The reciprocal is also immediate. 
4.3 Singularity of type III: higher order vanishing of A£ 

In this section we want to clarify the character of the phase transition in the case of a type III 
singularity. To keep things simple, let us restrict our attention to the quartic case (m = 2 in 
(34)), but the conclusions readily follow for the general situation. Assume that for t = T > 
the type III singularity occurs, without loss of generality, at the right endpoint of the support. 
This means that for t G (T — s,T + s), t ^ T, 

A(x) = (x — a\)(x — a 2 ), B(x) = (x — b)(x — b), 
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where a\ < a 2 and Imft > 0, in such a way that both Imfc and (Re b — a 2 ) are small with 
respect to 02 — cti, as t — > T. By (42), for < \t — T| < e, 



1 



b 



(ai - a 2 )B(ai)' ° 2 (a 2 - a 1 )B(a 2 ) ' 

1 ^ = ^ = _ 2 

(6 - b)A{b) ' (6 - 6) A(6) ' 



(92) 
(93) 



For the sake of brevity we denote 

a = ai, d = a 2 — a%, 5 = Keb — a 2 , v = (Im6) 2 > 0. 

Then 

A(b) = 5{d + 5)-v + i{25 + d)(Imb), \A(b)\ 2 = (5 2 + v) ({d + 5) 2 + v) 
B{ai) = {d + 5) 2 + v, B(a 2 ) = 5 2 + v, B( ai )B(a 2 ) = \A(b)\ 2 . 



(94) 



By (93), 



(Re b) 



Im A(b) 



2(Im6)|A(6)| I 



d + 25 



(Im6) 



Re A{b) 



v - 5(d + 5) 

2(Im6)|A(6)| 2 ~ 2(Im6)|A(6)| 2 



Combining it with (92), we arrive at 



d = d 2 — d\ 



2(5 + d) 2 + 5 2 + 2v 



d 



5 = (Re b) — d 2 = 

. _ v-5(d + 5) 
V ~ \Mb)\ 2 ■ 



\Mb)\ 2 

l(d + 5) 2 + ld(d + 5) + v 

d \aW 2 : 



(95) 
(96) 
(97) 



Observe that by (96), 5 < for t in a small neighborhood of t = T, so that we can use 5 as 
a new variable. 



Dividing ( 95 ) and ( 97 ) by ( 96 ) we get 

dd _ 
86 ~ 
Ov 



(5 + d) 2 + 5 2 + 2v 
\d + 5) 2 + \d{d + 5) + v'- 



8{d + 5)-v 



85 



\d + 5) 2 + \d(d + 5)+v' 



(98) 
(99) 



Notice that the right hand sides are analytic at 5 = and v = as long as d > 0. Thus, 
for any do > there exists an analytic solution (d(5), v(5)) of (98)-(99) satisfying the initial 
conditions 

d(0) = d , v(0) = 0. 
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The series expansion of this solution at 5 = gives 



d = d n - -<5 - — <r ^S 6 



3 9d 27d? Q 81dg 



3 9d 18d§ v ; 



Replacing it in ( 96 ) yields 

d_ 

dt 



(S 3 ) = 35 2 5 



27 



7 



d V 6do 



If we denote <?(i) = <5 3 (i), we will obtain from here that 



l-^ 1/3 (i) + 0(g 2/3 (i))) dq 



27 



or 



We can rewrite it as 



(t)-—q 4 / 3 (t) + 0(q 5 / 3 (t)) 



7 5 2 + 0{5 3 \ 



27 

8d~o 



dt, 



(t-T). 



and finally, 



24d 
7 



27 



(t-T), 



24d, 



o 



<5 Z + 0(5 6 



2d, 



1/3 



(t-T) 



1/3 



(100) 



(101) 



The analytic function of 5 in the left hand side of ( 101 ) is invertible, and straightforward 
computations show that 



S(t) 



2d, 



3 m (t - T) 1 / 3 + -^(t - T) 2 / 3 + 0(t — T), t -> T. 



32(£ 



(102) 



On the other hand, by (100), 



5^0, 



so that ( 102 ) yields 



Tmfc^ _^|t_ T | 1 /3 + _?L_( t _ T )2/3 + ( t _ T ) i t _> r . 



2d, 



1/3' 



32d, 



5/3 



(103) 



Formulas (102) and (103) show that both Re(6) — 0,2 and Im(6) have the asymptotic order 
0(\t — T\ 1 / 3 ) in a neighborhood of the critical value t = T at which b impacts ai and returns 
to the complex plane, i.e. where the type III phase transition occurs. This is consistent with 
the result obtained previously in [HJ Lemma 8.1]. 
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It is convenient to point out also a certain universal behavior this transition exhibits, 



expressed in the fact that the leading terms in ( 102 ) and ( 103 ) depend on the initial value 



da only. In particular, we see that b impacts ci2 from the complex plane with the asymptotic 
slope of 7r/6, both for the incidence and the reflexion angles, as it were for the actual reflection 
law in Optics. 

We switch now to the analysis of the Robin constant p(St); in the case we are studying, 
St = [0,1,0,2] in a small neighborhood of t = T, so that with d = a-i — a\, 

p(S t ) = 2log2-logd. 



From (95) we obtain that 
d 



dt 



p(S t 



d 
d 



2 (6 + d) 2 + 5 2 + 2v 



d 2 



\A(b)\< 



for < \t-T\ < e. 



Using that 5, v 



and 



as t — > T and expressions (|94j) and (100), we get that 
(5 + d) 2 + 5 2 + 2v 



lim 



d 2 



1, 



lim 

t-»T 



8 2 



A, 



d 



Taking into account (102) we conclude that 

2 



d_ 
dt 



p{S t 



t - T|~ 2 /3 (x + o(l)) , t->T. 



3 d. 



4/3 



(104) 



In other words, in the case of a singularity of type III, we have again a third order phase 
transition with an infinite algebraic jump of the third derivative of the free energy at the 
critical time, with the exponent —2/3. 



4.4 Birth of new local extrema 

We finally turn our attention to the situation created when a pair of complex conjugate zeros 
b and b of B collide at 60 G M \ St and become two new simple zeros 6_, b + of B. It was 
mentioned that all a^'s are analytic through t = T, and thus by (39), B and Rt are also 
analytic functions of the parameter t. 

Using the notation introduced above, our assumptions can be written as follows: for a 
small e > 0, there exist polynomials A, B and h, continuous with respect to the parameter 
t € (T — e, T + e), but represent, generally speaking, different real-analytic functions of t for 
t < T and t > T, such that 



A(x;t) = A(x;t), te(T-e,T + e), 

B(x . t) = H x - b )( x - l ) B & t )> t€(T-e,T), 

\(x-b.)(x-b+)B(x;t), te[T,T + e,T), 

h(x;t) = h(x;t), t E (T - e,T + e), 
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where b and b± are continuous functions of t such that 

b(t = T-) = b-(t = T+) = b+(t = T+) = b eR\ S t . 

Notice that b G C \ R for t G (T — £,T); without loss of generality, Im6 > 0. A priori, we 
do not assume that b± are real- valued for t G (T,T + e), so the two possibilities are either 
6_ < 6+ or 6_ = 6+ G C \ M. 

We denote = h(x)/(A(x)B(x)), omitting from the notation when possible the ex- 
plicit dependence on t. 



From ( 42 ) it follows that 

6-6 



i G (T — e, T), (105) 



5^' = ^ + doe) 



Subtracting/adding both equations in (105) we easily get that 

d //T ,, 9N _^ .,. d ^ 7 Imo(6) 

It follows in particular that 

6(t) = Vff(M(r-t) 1/2 (i + (i)), t^T-. 

Observe that these formulas show that the collision of b and 6 on R \ St can occur only at a 
position 6o where 

3(6 ) > 0. 

Moreover, assuming that for t > T the new zeros b± are complex conjugate, the same formulas 
apply. This yields the partial conclusion: the scenario when the complex zeros of B collide 
atM\St and bounce back to the complex plane can occur only when 60 is either a zero or a 
pole of g. This is the situation, for instance, when 60 coincides with one of the endpoints of 
St, and in this case we get a type III phase transition. 
Thus, let us assume that 

g(b ) > 0. (107) 

This is always the case, for instance, if St is a single interval. Then, 6_ < 6 + for t G (T, T+e). 
Denoting again 5 = (6+ — 6_)/2, we obtain in the same fashion that 

5(t) = N /fl(6^)(t-T) 1 /2(i + (i)) j t ^T + . 

Finally, we have mentioned that a singularity of type III is a limit case of the situation 
analyzed here, when bo coincides with one of the a^'s. Furthermore, under assumption (107), 
a collision of a pair of complex conjugate zeros of B at 60 G R\ St is followed by a type I phase 
transition (birth of a new cut). However, these two phenomena cannot occur simultaneously: 
if 6 6 (a2fc,a 2 fc+i), conditions 

&0 r a 2k+l 



a 2k J a 2k 



^Rdy)dy= / ^/R^y)dy = 



(see (24)) are incompatible. 
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5 The quart ic external field 



In this section we consider in detail a particularly important case of a quartic potential, 
i.e. when m = 2 in the representation ( |32[ ). Observe that this is the first non-trivial situation, 
since for m = 1 (quadratic polynomial) all calculations are rather straightforward. According 
to Section [3j for t > 0, St has the form (15) with either p = 1 ("one-cut regime" or "one-cut 



case") or p = 2 ("two-cut case"). Additionally to the description of all possible scenarios 
for the evolution of St as t travels the positive semi axis, we characterize here the quartic 
potentials <p for which St is connected for all values of t > 0, as well as those for which the 
singularity of type III (higher order vanishing of the density of the equilibrium measure At) 
or the birth of new local extrema occur. 

Roughly speaking, the evolution of St can be described qualitatively as follows: for the 
quartic potential ip there exists a two-sided infinite sector on the plane, centered at a global 
minimum of p and symmetric with respect to the horizontal line passing through this min- 
imum and with the slope 0.27872057. . . , such that St is a single interval for all values of 
t > if and only if the other critical points of p lie outside of this sector. Otherwise, the 
positive i-semi axis splits into two finite subintervals and an infinite ray. The finite subin- 
terval containing t = (which may degenerate to a single point t = 0) corresponds to the 
one-cut situation, and for the neighboring finite interval the support St has two connected 
components. Finally, the infinite ray corresponds again to the one-cut case. 



Recall (see Sections 3.3 and HI) that the transition from one to two cuts occurs always by 



saturation of the inequality in <\8h, while the transition from two to one cut occurs by collision 



of some zeros of the right hand side of ( 34 ) . 

Let us give the rigorous statements. For any quartic real polynomial <p with positive 
leading coefficient we define the value sl{(p) as follows: let (0 £ R denote a point where (p 
attains its global minimum on M (which can be unique or not), and £l G C, any other critical 
point of <p (zero of ip'). Then 



sl((p) 



ImCi 



, ifReCi/Co, 

Co - Re Ci 
+00, if ReCi = Co- 



Geometrically, sl{p) is the square of the slope of the straight line joining Co an d Ci- Notice 
that a real cubic polynomial has either 3 real zeros, or one real and two complex conjugate 
zeros, so that there is no ambiguity in the definition of sl(-). 

Next, we define the critical slope: let s = 0.077685 . . . denote the only positive root of 
the equation 

32s 3 - 17s 2 + 14s - 1 = 0. (108) 

Explicitly, 



96 

Alternatively, 



1 / r _ 1055 + \ w . th r= 3/ 5 / 3072> /g_3 107 \ >Q> 



10 -7 2 
30 ' 



(109) 
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where 

--II - - 24V5 ) " 3 + ( 59 - 24 ^) - 1 | = 2.7 6938 ... 




5 / 

is the only real solution of the equation 

4 7 3 + 157 2 - 200 = 0. 

Theorem 5.1. Let the quartic potential p be given, and s denote the critical value as described 
above. Then, 

Case sl(p) > s: St is a single interval for all values oft>0, no singularities occur. More- 
over, non-real zeros of Rt in (34) move monotonically away from the real line if and 
only if si ((p) > 1. 

Case sl((p) = s: St is a single interval for all values oft > 0, but there exists a (unique) value 
oft for which a type III singularity occurs (and this is the unique phase transition). 

Case < sl(ip) < s: St evolves from one cut to two cuts, and then back to one cut, presenting 
once the birth of new local extrema, a singularity of type I and a singularity of type II, 
in this order. No other singularities occur. 

Case sl{<p) = 0: if (p attains its global minimum at a single point, then St evolves from one 
cut to two cuts, and then back to one cut, presenting once a singularity of type I and a 
singularity of type II, in this order. No other singularities occur. 

If ip attains its global minimum at two different points, then St evolves from two cuts 
to one cut, and only a singularity of type II is present. 

We obviously consider t > 0; the value t = is not regarded as a singularity. 

Remark 5.1. Observe that for ip with more than one local extrema on M. there are no type 
III phase transitions. 

We can also easily characterize the quartic external fields p for which singularity of type 
III occurs (i.e. such that the zeros of <p' lie on the critical line) directly in terms of their 
coefficients. Indeed, if p>'(x) = x 3 + d2 x 2 + d\ x + do, then 

if d 2 \ 3 / d%\ 2d\ d±d 2 



So, without loss of generality, we may assume that 

ip'(x) = x 3 + dix + d . (112) 



Let xq be a real zero of tp'; according to Theorem 5.1 At will have a type III singularity for 
a certain value of t if and only if the other two zeros of <p' are of the form u± iv, u,v £ M, 
and the value 

v 

U — Xq 
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is a root of the polynomial (108). In particular, 



ip'(x) = (x - Xq) ((x - U) 2 + V 2 ) = (X - Xq) ((x - U) 2 + S(u - X ) 2 ) . 

Comparing ( 112[ ) and ( 113| ) we conclude that 

2u + xq = 0, u 2 + s(xq — u) 2 + 2xqu = d\, xo(d\ — 2xqu) = do. 



Eliminating xo and u from (114) we find that 

729dl s 3 - 81(4df + 9dl)s 2 + 9(-8d? + 27dg)s - Ad\ - 27d^ = 0. 



(113) 
(114) 
(115) 



The resultant of the polynomials in the left hand side of ( 108 ) and ( 115 ) is an integer multiple 
of (128c? 3 + 135o!q) . Since both polynomials share a common root, s, the resultant must 
vanish, and we conclude the following: 



Corollary 5.1. For an external field ip' with derivative of the form (112) the equilibrium 
measure Xt develops a type III singularity if and only if 

128dl + 135dg = 0. 



Using the substitution (111) we can easily extend this result to the general case: for an 
external field tp' such that 

<p'(x) = x 3 + c?2 x 2 + d\ x + do, 
the equilibrium measure Xt develops a type III singularity if and only if 

128df + 135dg - d 2 {d 2 {2d 2 2 - 9di) 2 + 2(16d?d 2 + 45d di - lOdodjj)) = °- 

For instance, direct substitution shows that (p'(x) = x 3 + 4x 2 + 2x — 8 satisfies this condition. 

According to ( 38 ) , in this case we will have a singularity of Type III for a finite value of 
t = T, where the density of the equilibrium measure vanishes with the exponent 5/2. The 
value of T can be found by the procedure described at the end of this section. 



Theorem 5.1 is a consequence of Theorems 5.2 and 5.3 below, where some additional finer 
results on the dynamics of the equilibrium measure as a function of t are established. 

Since the problem is basically invariant under homotopy, horizontal and vertical shifts in 
the potential ip, as well as mirror transformation x i— > —x of the variable, in the rest of this 
section without loss of generality we assume that 



(p(0) = = min{(£>(x) : x G E}, and p'(x) = x(x — a)(x — /3), 
with both a and (3 in the closed right half plane. We have 

ip(x) = - x 4 + t 3 x 3 + t 2 x 2 , t 3 



•-(a + /3), t 2 = -a/3, 



(116) 



(117) 



so that we may suppose that one of the following two generic situations takes place: 
Case 1: (p' has three real roots, and < a < (3 < 2a; 
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Case 2: ip 1 has one real root, at x = 0, a is in the first quadrant, and a = (3. 



Condition in Case 1 is equivalent to saying that ip has on M two local minima, at x = and 
x = (3, and a maximum at x = a, in such a way that = <p(0) < < (obviously, 

any general situation when ipt has three real roots can be reduced to Case 1 by an affine 
change of variables and by adding a constant to p). Since the situation (3 = 2a (or (p(/3) = 0) 
is equivalent to an even external field, in what follows we assume that the inequalities are 
strict: < a < /3 < 2a. 

Case 2 means that tp has only one local extremum on the real line. 

In this way, the only case excluded from the analysis is <p(x) = x 4 /4, for which sl(cp) = +oo 
and the situation is trivial. 



For t > the polynomial Rt in the right hand side of (34) has degree 6, and the identity 



(35) takes the form 



Rt(x) = (ip'(x)) 2 — 2t x 2 — dt x — et 



A(x)B 2 (x), (118) 
for certain constants dt,et £ K. The following technical result will be useful in what follows: 



Lemma 5.1. Assume that for Rt in (118), 



or equivalently, Rt(z) 
identities hold: 



A(z) = {z- a)(z - 
{z-a){z-b) 4 (z 



B{z) = (z-bf 



(119) 



c). Then if cp is given by (117), the following 



a + c + 4b = 2(a + (3), 
6b 2 + 46(a + c) + ac={a + /3) 2 + 2a/3, 
b(2b 2 + 36(a + c) + 2ac) = a/3(a + p), 
b 2 (b 2 + 46(a + c) + 6ac) = a 2 f3 2 - 2t. 



(120) 



Proof. This is a straightforward consequence of replacing (119) in (118) and equating the 
coefficients in both sides. □ 



Remark 5.2. It is important to observe that the identities (120) remain valid under a homo- 
thetic transformation 

a i->- qa, b h-» qb, c4 qc, q4 qa, (3 h-» q(3, t \-t q 4 t, 

for q > 0. In other words, a linear scaling in space yields a quartic scaling in time (or 
temperature). 

Let us consider first Case 1 (with strict inequalities). The theorem below is the quanti- 
tative description of the following evolution: for small values of temperature t a single cut is 
born in a neighborhood of the origin. The other two (double) zeros of Rt are real and close 
to x = a and x = (3, moving in opposite directions. At the first critical temperature t = T\ 
a bifurcation of type I occurs: the rightmost double zero splits into two simple real zeros, 
giving birth to a second cut in the spectrum. This configuration is preserved until both cuts 
merge at a quartic point at a temperature t = T2 (phase transition of type II). After that, 
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two complex conjugate double roots of Rt drift away to infinity in the complex plane, so that 
for the remaining situation we are back in the one-cut case. 

For the sake of convenience, here we use all introduced notations interchangeably, 

■ d A 

a = —a = Ooa. 

at 



Theorem 5.2. Let ip be an external field given by (116) with < a < (3 < 2a. Then there 
exist two critical values < T\ < T 2 such that: 

• (Phase 1 ): for < t < T±, A and B in ( 118| ) have the form 

A(x) = (x-ai)(x-a 2 ), B(x) = (x - b x )(x - b 2 ), (121) 
with a\ < < a 2 < b\ < b 2 , and St = [(21,02] (one-cut case). 

Parameters a k , b k are functions oft and satisfy the system of differential equations 

2 ■ 1 

A'(a k )B(a k ) A[b k )B'{b k ) 

with the initial values 

a x (t = 0) = a 2 (t = 0) = 0, b 1 (t = 0) = a, b 2 (t = 0) = /3. 
For0<t< Ti, 

di < 0, d 2 > 0, b\ < 0, 6 2 > 0, (123) 

and function 

F(t)= / y/Ajs)B(s)ds (124) 
Ja 2 (t) 

is monotonically decreasing and positive in (0,Ti), with 

F(0) = <p{f3) = (2 " ~®P 3 > and F(T\) = 0. 

• (1st Transition, phase transition of type I): for t = T\ such that F(T\) = we 
have S Tl \ S Tl = with ■= b 2 (t = Ti), which is a singular point of type I 
("birth of a cut"). 

• (Phase 2): for T\ < t < T 2 , A and B in ( |118[ ) have the form 

A(x) = (x — a±)(x — a 2 )(x — a^)(x — 04), B(x) = x — 61, (125) 

with a\ < < a 2 < b\ < a% < b^ 1 ' < 04, and St = [a±, a 2 ] U [03, 04] (two-cut case). 
Moreover, these values satisfy the system of differential equations 

a k = 777 r , fc = l,...,4, and 6! = ——. 126 
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Parameters ai(t) and 02 (i) are continuous att = T\, while 

a 3 (t = Ti) = a 4 (t = Ti) = £(t = Ti) = b«. 
TTie uaZue of ^ € (02,03) *s determined by 

r ° 3 x-C 



(127) 



«2 



dx = 0. 



In particular, 

01,03 < 0, and 02,04 > 0. 

T/ie critical temperature T2 is determined by the collision condition 0,2 (t = T2) = 03 (i 
T 2 ). 



f^nrf Transition, phase transition of type II): for t = T2, A and B in (118) have 
the form 

A(x) = (x- ai)(x - a 4 ), B(x) = {x - b (2) ) 2 , 



where 



= lim 02 = lim 03 = lim ( = lim b\, 

t^T 2 - *->T 2 - t->T 2 - t^T 2 - 



so that ai < < b^ < b^ < 04, and St = [01, 04]. Points a±, b^ and 04 can be found 
from the first three equations in (120) (with a = a\, b = b^ and c = 04), while the 
value 0/T2 is obtained from the fourth equation (120). 

Since b& G S T2 , x = b^ is a singular point of type II, and Rt 2 has a root of order 4 
atz = b&>. 



(Phase 3): for t > T2, A and B in (118) have the form 

A(x) = {x-a 1 )(x-a 4 ), B(x) = (x - b x )(x - b{), (128) 

with a± < < b^ < b^ < 04, Im&i > 0, and St = [01,04] (one-cut case). 

Moreover, these values satisfy the system of differential equations of the form (122) 
(setting now 62 = b\); in particular, \m.b\ grows monotonically with t, and 



lim Imbi = +00, lim Rebi 

t— >+oo t— >+oo 



lim 

t— >+oo 



ai + 04 a + f3 



(129) 



Remark 5.3. If (3 = 2a, then F(0) = 0. The evolution is as described, except for T\ = 
(Phase 1 is missing). 



Proof. Using (40) and representation (27) we conclude that there exists e > such that for 
< t < e we are in the one-cut case, and formulas ( 121 ) hold, with both b\ and 62 close to a 
and (3, respectively. In this situation, h = 1, so that (122) is a particularization of (42), and 
the inequalities (123) are just straightforward consequences of (122). Observe also that (24) 
implies that 



l>2 



y/A(s)B(s)ds > 0. 



(I 2 
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Furthermore, for any finite initial positions a&(i = 0) := a?, bk(t = 0) := b®, with 
< < b® < 62, the solution of the system of differential equations (122) exhibits collision 
(of a 2 and bi) in finite time. Indeed, by (123), for t > (and before the collision), 

a± < a\ < a!] < a 2 < b\ < b\ < 6° < b 2 . 

Since 62 < (b 2 — &?) _3 > solving the corresponding ODE we conclude that 

b 2 <G 1 (t) := & o+(4t + (^-6?) 4 ) 1/4 . 

Analogously, 



01 > G 2 (i) := a§ - (8i + (a§ 
Replacing these bounds in ( 122| we get 

2 

a 2 > 



or 



(61 - a 2 )(6'i 



b, < 



0-2) 



(b° 1 -G 2 (t))(G 1 (t)-a° 2 )(b 1 -a 2 ) 7 
1 

(6?-G 2 (i))(G 1 (t)-a0)(a 2 -6 1 )' 



1 

2 eft 



(61 - a 2 ) 2 < 



(6?-G a (t))(Gi(t)-aS)" 



Thus, a collision will occur by time t = T if 



o (6S-G2(0)(G!(t) 



(6? 



,0\2 



6 



Since Gi(t) ~ t 1//4 , G 2 (t) ~ t 1 / 4 as t — > 00, the integral in the left-hand side diverges as 
T —> +00, which proves that there will always be a collision in a finite time T. Observe that 
function F in (124) is well-defined in the whole interval (0, T), and that the integrand in 
(124) is, up to a constant, the analytic continuation of the density of At, see (38). Using (26) 

la 2 {t) \fMsj 



we conclude that 



F'(t) 



< 0. 



However, at the collision time T, 

-6 2 



F(T) 



y/(x — a\){x — a 2 )(x — b\){x — b 2 ) dx < 0, 



a 2 



which shows that there is a unique time T\ < T for which F{T\) = 0. 

From the positivity of the measure At and expression (36) it is easy to conclude that all 
roots of Rt, which are not endpoints of the support, need to be double. Hence, for t > T\, 
the rightmost double root b 2 splits into a pair of simple real roots a% and 04, giving rise to 

x 



formula (125). We apply again Theorem 3.1 with h{x) 



£, which yields (126). 
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Observe that it is not straightforward to deduce the sign of b\ from 



Taking into account the initial values (127) we see that for a small e > and T\ < t <T\+e, 
Q is close to CI3, so that b\ < Q. This implies that immediately after the birth of a cut (t = T±), 
b\ < 0, so that point b\ still moves to the left "by inertia". In that range of time, 02 and 
03 (and in consequence, also b\ and C) are in the collision course, and collision occurs in 
finite time, when all these four points merge simultaneously. This critical time T2 > T\ can 
be characterized by the appearance of a quadruple root b^ of Rt 2 inside St 2 , so that the 



system (120) is valid 



Taking into account the monotonicity of St , we see that for t > T2 the quadruple root b^ 



of Rt splits into two complex-conjugate roots 61 and b\, and formulas (128) hold. 

hi) 



Adding the equations (42) for b^s (62 



d 
df 



(Im&i) 2 



Re 



we obtain 
1 ReA(h) 



(130) 



A(h) \AQn)\* 

Observe that Re^4(2;) = is an equation of an "East- West opening" rectangular (or equilat- 
eral) hyperbola r = Tt with its vertices at a% and 04, and {z G C : HeA(z) < 0} corresponds 
to the connected component Ut of its complement in C containing the segment joining a\ and 
04. 



r*\ 






ai 1 




04 







Figure 1: Hyperbola Tt and the domain Ut (shaded). 



Hence, we conclude from (|130) that 

Im(6i) > ^ b% G Ut- 



d_ 
dt 



(131) 



Moreover, the monotonicity of the support St (or equivalently, the fact that d\ < < a'4) 
implies that 

tx < t 2 => U n C U T2 . 
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Clearly, for t > T2, h\ € Ut, so that in this range of t, Im(£>i) grows monotonically, and 
we have a one-cut case for all t > T^. 

Finally, since by assumptions (116), (p(x) > for all x € M\ {0}, using the representation 
(27) we conclude that 



lim ( 

t-H-oo 



-01 ] 



lim 04 = +00, I St 



t—t+OO 



t>0 



get 



Regarding the second limit in (129), observe that from equations (42), taking b 2 = b\, we 
d ( a\ + 04 



2 Re 61 — (ai + 04) 



dt V 2 ) dt y ±! \Mh)\ 2 

An immediate consequence of these identities is that the centers of masses of the zeros of A 
and of the zeros of B are always in a collision course. 

A comparison of the coefficients at x 3 ,x 4 and x 5 in both sides of (118) yields the system 

4 Re 61 + ai + a 4 = 2 (a + /3), 
6(Re 61) 2 + 2(Im &1) 2 + 4(ai + a 4 ) Re h + a x a A = (a + /3) 2 + 2a/3, (132) 
4|6i| 2 Re6i + 6 ((Re6i) 2 + 2(Im6i) 2 ) (ai + a 4 ) - 4aia 4 Re 61 = -2a/3(a + /3). 



An assumption that a\ + 04 = 2 Re 61 (collision) in the last two equations in ( 132 ) implies 
that 2a 2 + 2/3 2 — 5a/3 = 0, which is possible only if /3 = 2a. This is the symmetric case not 
considered here. 

Hence, we conclude that in our situation there exist the limits 



lim 

t— >+oo 



d\ + a 4 



lim Re 61 . 

t—t+oo 



By using this in the first identity in (132) we conclude the proof of (129). 
Next, we turn to Case 2. 



□ 



Theorem 5.3. Consider the external field given by (117) with (3 = a, and with a in the first 
quadrant. Then: 



(a) if 



arg(a) > arg(l + iy/s), 



where y/s = 0.27872057. . . , and s is the only positive root of the equation (108), then 
for all t > 0, polynomials A and B in (118) have the form 

A(x) = (x-a 1 )(x-a 2 ), B(x) = (x - h){x - h), lm&i>0, (133) 

with a± < < a,2, and St consists of a single interval [01,02] (one-cut case). 

//arg(a) > 7r/4, then functions —a\, 02 and Im h\ grow monotonically with t from their 
initial positions 0, and Im a, respectively, to +00. Otherwise (i.e., i/arg(l + iyfs) < 
arg(a) < vr/4^, there exists a critical value Tq > such that Im&i is monotonically 
decreasing for < t < To, monotonically increasing for t > Tq, and Im(&i(t = To)) > 0. 
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(b) if 



< arg(a) < arg(l + iy/s), 
then there exist three critical values < To < T± < T2 such that: 



for < t < Tq, polynomials A and B in (118) have the form (133) with ai < 



< 02; Imfci is monotonically decreasing for < t < Tq. At the critical value Tq, 
f)W := 5 1 (^ = 7q) g R \ (01,02), which is a birth of new local extrema; without 
loss of generality, we assume > ci2(t = Tq). During this phase, St consists of 
a single interval [01 , 02] ■ 



for Tq < t < T\, polynomials A and B in (118) have the form (121), with a\ < 
< 02 < b\ < &2; and St = [01,02] (still one-cut case). 

At t = T\, function F defined in (124) satisfies F(Ti) = 0, and 62 = T\) = 



is a singular point of type I ("birth of a cut"). 

forT\ < t < T2, Phase 2 described in Theorem \5.2\ takes place. Second transition, 
also described in Theorem 5.2, occurs at t = T2 by the same mechanism (collision 



of four zeros of Rt). This situation corresponds to the two-cut case. 

finally, for t > T2 we have Phase 3 as described in Theorem \5.2\ We are back in 
the one-cut case. 



Proof. As in the proof of Theorem 5.2 using (40) and representation (27) we conclude that 



there exists e > such that for < t < e we are in the one-cut case, and formulas ( 128 ) hold. 



Using the equivalence ( 131 ) we conclude, in particular, that if a € Uq, that is, if | Re a\ < Im a, 
then Im(6i) grows monotonically with t, varying from Im(a) to +00 as t travels (0,+oo). 
Clearly, in this situation we have a one-cut case for all t > 0, which establishes the first part 
of (a). 



As it was discussed before, for a fixed A, the solutions bi and b\ of (130) exhibit collision 



in a finite time. The second part of Case (a) corresponds to the situation when for a To > 
before the collision, b\ 6 r^ (see Figure [l]). 

Finally, consider the case when 61 hits the real line in a time To. Clearly, for t < Tq, 
b 1 <£ U t , so that M°) = bi(t = Tq) G K \ (a$,aP,), where a° k := a k (t = T ). Assume first that 
b(°) ^ a ° k , k = 1, 2. At t = Tq we have that A and B in (fUSl) have the form 



A(x) = (x — a\)(x — 02), B(x) = (x — &1) 2 , 

with (without loss of generality) a± < < 02 < b\. As t becomes greater than To, the only 
possibility is that b± splits into two real zeros, so we are left in the situation of Phase 1 of 



Case 1 (see Theorem 5.2). From this point the evolution of A, B and St follows exactly the 
Case 1. 

The boundary between (a) and (b) is precisely when 

&i(t = T ) = a 2 (t = T ), 



that is, when Rt has a zero at b\ of order 5; using the expression in (118) we get that for 

t ~ 'lh, 

Rt(x) = x 2 (x — a) 2 (x — a) 2 — 2t x 2 — dt x — et = (x — oi)(x — 61) 5 . 
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We can use the identities ( 120 ) from Lemma 5.1 with a = a%, b = c = b\ and (3 = a. Moreover, 
by the homogeneity of these identities (see Remark 5.2) we may assume that a = 1 + i\/s, 



with s > 0: the boundary between case (a) and case (b) will be given by the ray ua, with 
u > 0. 



Thus, from (120) we obtain 



ai + 56i = 4, 
56i(ai + 2&i) = 2(3 + s) 
56?(ai + h) = 2(1 + s) 



2t 



:i + s) 



56?(2ai 



These equations must have a solution with a\ < < b\ for t > 0. Solving the first two for ai 
and b% and replacing the solutions in the third equation (preferably, using a symbolic algebra 
software) we obtain 

200 



4r 3 + 15r 2 



0. 



y/W - 30s 2 , 



which is equivalent to (108). Direct calculations show that the discriminant of the cubic 
polynomial above is negative, so it has a unique (positive) root s = 0.077685 . . . , which 



satisfies identities (109)-(110), and two complex conjugate roots. Finally, a substitution of 
the obtained values for s, a%, b\ in the last equation renders t = 0.0339206 



> 0. 



□ 



We conclude this section with a number of remarks that we consider relevant. 



The results from Theorem 3.3 allow us to explore the dependence of the support St (and 
of the corresponding equilibrium measure At) from the rest of the parameters of the problem, 
in our setting, from the position of the critical points a and (3. For instance, when St consists 
of a single interval ("one-cut case") differential relations (54)-(55) apply. Observe e.g. that 



in the case when the center of masses of a and (3 increases (hence, t% in (117) decreases), 
both endpoints of St move to the right. 

In a certain sense, the considerations above show that the most general situation corre- 
sponds to Case 2 with a lying in the complex plane but below the critical ray emanating 
from the origin and passing through 1 + iyfs (or, in the terminology introduced at the begin- 
ning of this section, when < sl((p) < s). This and only this situation exhibits all possible 
transitions occurring in the quartic case. 

Observe that the external field ( 117| ) with (3 = a, and with a in the first quadrant is 
convex if and only if arg(a) > tt/6 > arg(l + i^/s), where s is the critical value given by 



(109)-(110). Hence, the persistence of the one-cut case for all t > under assumptions of 



convexity (see e.g. [72]) is a consequence of the statement above. 

Example 5.1. Let us return to the analysis of Bleher and Eynard in [S], where they consid- 
ered the external field of the form 



x 

T 



4ciar 



+ (2c? - l)x 2 + 8cix, with ci G (-1,1) 



(134) 



This family always exhibits a second transition (see Theorem 5.2) or a singular point of 
type II for t = T 2 = 1 + 4c?, in such a way that -2 = at{t = T 2 ) < 6 (2) < a 2 (t = T 2 ) = 2. 
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According to our analysis, for c\ ^ there is also a first transition (birth of a cut) at T\ < T2, 
not described in [9], but predicted in [36J and found numerically in pQ. 

Expression ( |134 ) is an alternative parametrization of the quartic external field; the range 
ci E (— 1, 1) covers all the cases when a two-cut St arises, described in Theorem 5.1 However, 



with the purpose of characterizing the existence of phase transitions we find the description 
given in our Theorem |5.1| more transparent. 

The analysis in [9] was extended in |55j to the case when tp is a polynomial of even degree 
(greater than 4) such that in the critical time t = T2 the density of the equilibrium measure 
has a zero of higher order 2m, m > 1. 



We finish this section observing that using the system of equations (120) we can find 



algebraic equations for the quadruple critical points and and for the corresponding 
endpoints of the support S T ° and S T2 , respectively (so that these values are algebraic functions 
of a and /3, and hence, of the coefficients of the external field). Indeed, solving the first two 



equations in (120) for a and c, we get 

a + P 
a + P 



a 

c 



2b - y/4b(a + P)- 2aP - 6b 2 , 
2b + y/4b(a + P)- 2a(3 - 6b 2 , 

and replacing it in the third one yields 

p(b) = Wb 3 - 10(a + P)b 2 + 2((q + P) 2 + 2aP)b - aP{a + P) = 0. 

If < a < P < 2a, this is equivalent to finding real roots of the polynomial 



P [X) 



10x J 



10(1 + u)x 2 + 2((1 + u) 2 + 2u)x - u(l + u) 



(135) 



(136) 



(137) 



for the values of the parameter u £ (1)2]. Its discriminant (easily found with the help of a 
computer algebra system), 



20 (4u 6 



12n 5 + 29n 4 



38u 3 + 29u 2 - I2u + 4) , 



is positive for u G (1,2], and the roots of p* for u = 1 are (5 ± \/5)/10 and 1. Hence, this 
polynomial has three positive simple roots in the indicated range of u. However, only the 
middle one guarantees that 



2x{\ + u) 



u 



3x z > 0, 



which, according to (135), is a necessary condition for having real solutions a and c in (120). 



As observed, for u = 1 this root is (5 + \/5)/10 ~ 0.7236, while for u = 2 it is 1. 

Hence, in the case when < a < P < 2a, the recipe for finding the (unique) critical point 



is as follows: with u = P/a £ (1,2], find the root x* of (137) lying in 

[5 + V5 



10 



1 



and take M 2 ** = ax* . Then the endpoints of the support are obtained by replacing b = b^ in 
(135), and taking a±(t = T2) := a < b^ < a,2(t = T%) := c; the value of T2 is finally computed 
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from the fourth equation in (120). In this case, as we have seen, x = is a singular point 
of type II (zero of the density of the equilibrium measure). 

Let us consider the case when a is in the first quadrant, and a = (3. 
Im(a)/Re(a) > and dividing (p6| by (Re(a)) 3 , 



Denoting £ 
we arrive at the polynomial equation 



P {X) 



5x 6 



Wx 2 + 2(3 + C)x 



(i + C 2 



o. 



(138) 



whose discriminant, as natural, is given by the left hand side of (108) with s replaced by s 



Thus, for £ > \/s, with s defined in (109)— (110), polynomial p* has only one real root, which 



yields non-real values of a and c in (135). In this case, no quadruple critical point appears. 



If on the contrary, < £ < y/s, then p* has three positive roots; only two of them give 



real values of a and c in (135). 



Summarizing, in the case when < £ := Im(a)/Re(a) < y/s, the recipe for finding two 
critical points is as follows: find the real roots x* of (138), take b = Re(a) x* and replace it 

JO = b (0) and c 



in (135). If a < c < b, then a = a(t = T 



b(t 



of To is computed from the fourth equation in ( 120| ). 

If on the contrary, a < b < c, then a = a(t = T2), b = b(t = T2) - 
the value of T2 is computed again from the fourth equation in (120) 



c(t = To); the value 
b^ and c = c(t = T 2 ); 
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